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Monica A. De La Paz 



I, Jerry Bryant, do declare that: 

1 , I am a United States citizen residing at 6861 Staffordshire St. Houston. IX 77030. 



2. I currently hold the position of Chief Technology Officer, Mead of Scientific Evaluation* 
Division of Business Development, Cell>Point, LLC, 7120 E, Orchard Road, Suite 350, 
FnglewmxL CO 80111. A copy of my curriculum vitac, including a list of my 
publications, is attached as Appendix A. 

3, I tiTu skilled in the synthesis and use of radk>nuclide-labclcd imaging agents, as evidenced 
by the following: 

• T have worked as Chief Technology Officer of Cell>Point, LLC since 2001. 

• T have been employed as the Director of Research Development of VeriMcd Research 
Corporation in Houston, TX from 2002-2004, 
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• My duties at Cell>Point and VcriMed include participation in the development of novel 
radiopharmaceuticals lor different diseases, such as cancer, cardiovascular disease, and 
diabetes, i have also been involved with studies directed to the synthesis and use ol" 
radionuclides, and studies directed to understanding the mechanism and biochemistry of 
the agents as it pertains in the pharmacokinetics and biodtstribution of the agents in 
animals and humans. In addition, T have been involved in the evaluation of new 
technologies for the treatment of cancer, cardiovascular disease, and diabetes. 

• 1 was employed as Chief Scientific Officer of Allcurc. inc., in Houston, "IX from May. 
2001 -May, 2002. My duties at Allcurc included designing and marketing in vitro and in 
vivo services to drug and bio-lech companies for evaluating compounds of interest 
(nuclear medicine, SCID mouse in vivo service and mechanism studies). 

• I have experience as a Research Assistant TI in the Department of "Nuclear Medicine, the 
Division of Diagnostic Imaging, of the University of Texas M.D. Anderson Cancer 
Center, from October, 2002 to June, 2003. 

• i am a co-inventor of two patent applications that pertain to imaging technology, 
including USSN 10/703,405 C k F4hyle:nedicysteine (FC)-Drug Conjugates, Compositions, 
and Methods for Tissue Specific Disease Imaging," Yang el a//); and USSN 10/942.615 
(^Mechanism-Based Targeted Pancreatic Beta Cell Imaging and Therapy," Yang el ciL), 

• I have been involved in funded studies pertaining, to the development and evaluation of 
imaging agents during the past five years, including: (i) a study of CT and MR1 
functional agent development and evaluation supported by VcriMed Research 
Corporation; (2) a study of 99m-Tc-Hthylenedicystcine (HQ-Drug Conjugates for Tissue 
Specific Disease Imaging supported by Cell>Point T f JLC; and (3) a study to compare Tc- 
90n>EC.Vdeoxyglucose (EC-DG) and FDG-PET scans for the evaluation of patients 



suspected of having persistent/recurrent squamous cell carcinoma of the larynx after 
definitive treatment with radiation therapy and the evaluation of primary lung cancer 
patients, sponsored by Coll>PoinU IXC See page 7 of Appendix A, 

• 1 am a co-author of seven articles and numerous abstracts pertaining to the evaluation and 
testing of radiolabeled imagine agents. See pages 9-14 of Appendix A. 

• I am also a co-author of two book chapters pertaining to radiolabeled imaging agents and 
their uses in chemistry and nuclear medicine. See page 1 4 of Appendix A. 

• Regarding my formal education, I have a Master of Science degree (1991) in 
Microbiology & Cell Science and Molecular Biology from the University of Florida, and 
a B.S. degree (1987) in Chemistry and Biochemistry from Tennessee State University. 

• [ have extensive experience in cell and molecular biology, as delineated in my curriculum 
vitae. See pages 2-3, Appendix A. 

4. ! have read the above-referenced patent application and the present Office Action, and am 
familiar with the technology. I understand that the examiner has rejected the claims 
because the examiner is of the opinion that at around November, 2002 when the- 
application was filed, a person of ordinary skill in the field of the invention would have 
been motivated to substitute the t8 F-Iabeled penciclovir probe of Iyer ef aL (J. Nucl. 
Med., 2001. 42, p. 96-105) with 99m~Tc-ethylcncdieystdnu {F.C)-aminopenciclovir to 
practice the method of Iyer et aL to monitor expression of HSVl-thymidine kinase. The 
examiner argues that the teachings of Zareneyrizi et aL (Anti-Cancer Drugs, 1 999, 1 0. p. 
685-692) and Yang et ah (Ann. Nucl. Med. ScL 2000, 13, p, 19-36) provide motivation 
to convert the penciclovir to amino penciclovir, and to label the aminopencielovir with 
99m-Tc-FC 



A person or ordinary skill in the art would not have been motivated to substitute the lH l ; - 
labeled penciclovir probe of Iyer et al. with 99nvrc-KC-aiiiiiiopcricic.lovir because it is 
highly unlikely that 99m-l*c-EC-ammopencielovir wonld be suitable for imaging ITSV1- 
Thymidine kinase reporter gene expression since it would not he a substrate for HSV1- 
thymidine kinase. 

This is the case because the prior art teaches that the acyclic side chain of auanoside 
analogues such as penciclovir and acyclovir were known to substitute for sugar moieties. 
The side chain oFTc-KC-amim^penciclovir does not structurally resemble a sugar moiety 
or any part of a sugar moiety by virtue of the inclusion of the amino modification of 
penciclovir and by virtue nf binding of EC (a molecule which docs nut structurally 
resemble a sugar moiety or pari of a sugar moiety) to the amino moiety of 
aminopcnciclovir. Information concerning the importance of structural similarity of the 
side chain to sugar moieties can be found, for example, in FJion et aL Proc. Natl Acad, 
ScL USA, Vol 74. No. 12, pp, V/lf^T/ft), 1 077 (Inhibit 1), Schacflfar et aL, J. Med. 
Chcm, 14, 367-309, 1971 (Kxhibit 2), and Schaeffer et aL Nature, 1978 Apr 
1 3;272(5654);583-5 (Exhibit 3). 

Furthermore, a number of acyclic guanosine analogues were known m the field at the 
lime of the filing date of the present patent application. These structural analogues were 
known to have acyclic side chains that include at least a portion of a sugar moiety. A 
number of references provided examples of acyclic ^uanosine analogues that were known 
to be substrates for HSV-thymidine kinase thai have side chains resembling parts of sugar 
moieties, including De Clercq et aL Nucleosides Nucleotides Nucleic Acids, 2001 Apr- 
Jul;20(4-7):27l-85 (Exhibit 4), ilsley et aL Biochemistry. 1995 Feb 28;34(8):2504-10 
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(Exhibit 5), Golbraikh et aL Nucleic Acids Res. 1 989 Oct U ;1 7(1 9):7965-77 (Exhibit 5); 
and Martin et aL, J Med Chcm. 1 986 Aug;29(S): 1 384-9 (Exhibit 6). 

8, Given this information concerning acyclic guanosinc auaiogues thai way known in the 
field, a person of ordinary skill in the Held of the invention would not have been 
motivated to employ Tc-EC-aminopeneielovir in the method of Iyer el to probe for 
I-iSV 1 ^thymidine kinase expression. 

9. I heieby declare that all statements made herein of my knowledge are true and that all 
statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of 
the United States Code and that such willful false statements may jeopardize the validity 
of the application or any patent issued thereon. 
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CERTIFICATE OF ELECTRONIC TRANSMISSION 
37C.F.R. § 1.8 

I hereby certify that this correspondence is being 
electronically filed with the United States Patent and 
Trademark Office via EFS-Web on the date-below: 



Veb on the date-i^elow: 
Monica A. Do La Pa z 



I, Jerry Bryant, do declare that: 

1 . I am a United States citizen residing at 6861 Staffordshire St. Houston, TX 77030. 



I currently hold the position of Chief Technology Officer, Head of Scientific Evaluation, 
Division of Business Development, Cell>Point, LLC, 7120 E. Orchard Road, Suite 350, 
Englewood, CO 801 1 1. A copy of my curriculum vitae and information concerning my 
expertise in the synthesis and use of radionuclide-labeled imaging agents, including a list 
of my publications, has previously been set forth in paragraph 2 and Appendix A of my 
first declaration (filed with the Patent Office on August 28, 2009). Cell>Point is a 
licensee of the technology of the present patent application. 



3. I understand from the Office Action dated December 4, 2009 that the patent examiner is 
of the opinion that my previous declaration filed with the Patent Office on August 28, 
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2009 does not present sufficient proof that it is highly unlikely that m Tc-EC- 
aminopenciclovir would be suitable for imaging HSV-1 thymidine kinase reporter gene 
expression since it would not be a likely substrate for HSV-1. 

4. While I disagree with the Examiner's opinion that I have not provided sufficient proof to 
support my opinion, in an effort to provide clarification I will provide an even further 
detailed explanation to make it easier for the Examiner to understand why it is unlikely 
that 99m Tc-EC-aminopenciclovir would be suitable for imaging HSV-1 thymidine kinase 
reporter gene expression. 

5. By way of background, thymidine kinase catalyzes the transfer of a y-phosphoryl moiety 
from ATP to 2'deoxythymidine (dThd) to produce thymidine 5' -monophosphate (dTMP) 
(reviewed in Tung and Summers, Amtimicrobial Agents and Chemotherapy, Sep. 1994, 
p. 2175-2179 (Exhibit A) and Omari et al. 9 BMC Structural Biology 2006, 6:22 (Exhibit 
B)). The reaction scheme is depicted in Exhibit C. As can be seen, thymidine kinase 
catalyzes the incorporation of a phosphate moiety at the 5-position of the deoxyribose 
moiety of thymidine. 

6. The available literature at around the time of the filing date of the present patent 
application teaches that certain molecules that structurally resemble the natural substrate 
of thymidine kinase can undergo phosphorylation by thymidine kinase. Some of these 
molecules, such as penciclovir and acyclovir, include an acyclic side chain. The acyclic 
side chain includes hydroxyl moieties are positioned in a manner that resembles a portion 
of the deoxyribose moiety of thymidine and in particular that portion that includes the '5 
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hydroxyl moiety. FPCV is another such example. The hydroxyl moiety of the side chain 
of FPCV structurally resembles a portion of the deoxyribose moiety that includes a 5' 
hydroxyl moiety (see Exhibit D, which compares the acyclic moiety of FPCV to the 
deoxyribose moiety of thymidine). 

99m-Tc-EC-aminopenciclovir would not be expected by one of ordinary skill in the art to 
be suitable for imaging HSV1 -thymidine kinase reporter gene expression because the 
available literature suggests that it would not be a good substrate for HSV-1 thymidine 
kinase. De Winter and Herdewijn (J Med Chem. 1996 Nov 22;39(24):4727-37; Exhibit 
E) teaches that molecules other than the natural substrate thymidine that are substrates of 
HSV-1 thymidine kinase are structurally strikingly similar to thymidine. See page 4727 
left column and page 4733. 99m-Tc-EC-aminopenciclovir, with its bulky 99mTc-EC 
moiety and the amino group, is not structurally similar to thymidine. Further, De Winter 
and Herdewijn teaches that in designing high-affinity nucleoside substrates of HSV-1 
thymidine kinase, "care should be taken to maintain the geometry of the base moiety and 
sugar hydroxyl functionalities." See abstract. A chemical drawing of 99m -Tc-EC- 
aminopenciclovir is depicted in Exhibit F. In comparison to the chemical diagram of 
thymidine in Exhibit C, there is loss of a hydroxyl moiety relative to thymidine (and also 
relative to penciclovir). Thus, in 99m-Tc-EC-aminopenciclovir the geometry of sugar 
hydroxyl functionalities has clearly not been maintained. These findings are in contrast 
to FPCV, which does not include a bulky group attached to the acyclic moiety of FPCV 
and which maintains conformation of the hydroxyl group of the acyclic moiety in the 5' 
position. 
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8. In view of the above, a person of ordinary expertise in the synthesis and use of 
radionuclidie-labeled imaging agents would not have expected that 99m Tc-EC- 
aminopenciclovir would be a suitable substrate for HSV1 thymidine kinase or a suitable 
agent for imaging HSV-1 thymidine kinase reporter gene expression. 

9. I hereby declare that all statements made herein of my knowledge are true and that all 
statements made on information and belief are believed to be true; and further that these 
statements were made with the knowledge that willful false statements and the like so 
made are punishable by fine or imprisonment, or both, under Section 1001 of Title 18 of 
the United States Code and that such willful false statements may jeopardize the validity 
of the application or any patent issued thereon. 



Date Jerry Bryant 
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Exhibit A 



Antimicrobial Agents and Chemotherapy, Sept. 1994, p, 2175-2179 
0O66-4804/94/$04,00+0 

Copyright © 1994, American Society for Microbiology 



Vol. 38, No. 9 



Substrate Specificity of Epstein-Barr Virus Thymidine Kinase 

PETER P, TUNG 1 and WILLIAM C SUMMERS 1,2 * 

Departments of Therapeutic Radiology 1 and Molecular Biophysics and Biochemistry, and Genetics, 2 
Yale University School of Medicine, . New Haven, Connecticut 06520-8040 

Received 18 January 1994/Retumed for modification 25 March 1994/Accepted 9 June 1994 

Purified recombinant protein encoded by the BXLF-I open reading frame of the Epstein-Barr virus genome 
bas thymidine kinase activity. The substrate behaviors of various nucleosides toward this enzyme were tested. 
Halogenated deoxyuridiues, zidovudine, and bromovinyJdeoxyuridine are efficient substrates, while acyclovir 
and dihydroxypropylmethylguanine are relatively poor substrates for the Epstein-Barr virus thymidine kinase. 



Like most other herpesviruses, Epstein-Barr virus (EBV) 
encodes a virus-induced protein with thymidine kinase (TK) 
activity. This enzyme catalyzes the transfer of the gamma 
phosphoryl group of ATP to the 5' hydroxy! of a variety of 
deoxynucleosides to produce the corresponding nucleoside 
monophosphate. It is presumed that the further phosphoryla- 
tions to the diphosphate and triphosphate forms are carried 
out by cellular kinases. In the case of herpes simplex virus type 
1 (HSV-1) and HSV-2, the substrate specificity of the viral TK 
is much less restricted than that of the mammalian cell TK 
(12). Thus, in infected cells the viral enzyme initiates the 
conversion of nucleoside analogs to phosphorylated forms 
which either are inhibitors of the viral DNA polymerase or are 
relatively toxic when incorporated into the viral DNA (5). If 
the EBV TK has a similar broad substrate specificity, it may 
provide a route to the activation of nucleoside analogs for the 
production of useful antiviral nucleotides. To examine this 
possibility, we examined the use of various nucleoside analogs 
by EBV TK made in a bacterial expression system programed 
with the cloned EBV DNA sequence encoding this enzyme. 

Reagents, Nucleoside analogs, thymidine, ATP, raorpholino- 
ethanesulfonic acid (MES), and bovine serum albumin (BSA) 
were purchased from Sigma, ( 14 C]thymidine and [-y- 32 P]ATP 
were obtained from Aroersham. [ 3 H]bromodeoxyuridine 
([ 3 H)BrdU) was purchased from New England Nuclear, Water 
was from a Millipore MilliQ system. All other compounds were 
reagent grade or better. 

EBV TK. The DNA sequence from nucleotides 144862 to 
142745 of the B95-8 strain of EBV was cloned into the pET3a 
expression vector (20) at the Ndel and BamHl sites. This 
vector was introduced into Escherichia coli SY 211 (25), a 
TK-deficient derivative of E, coli BL21 (DE3) (20), so that the 
only TK activity in the cell comes from the expression of the 
EBV sequences. Cultures were grown at 22°C in 4 liters of LB 
medium containing 200 jxg of ampicillin per ml. Once the 
cultures reached the mid-log phase of growth, isopropylthio- 
galactoside was added to 0,4 mM to induce the expression of 
EBV TK. Induction was allowed to continue for 16 h at 22°C. 
Twenty grams (wet weight) of cells was harvested, washed with 
STE (20 mM NaCI, 1 mM EDTA, 10 mM Tris [pH 8.0]), and 
resuspended in 50 ml of 10% sucrose-50 mM Tris-HCl (pH 
7.5) before being frozen in liquid nitrogen, The frozen cells 
were defrosted in 100 ml of 10% sucrose-200 mM NaCI-30 



* Corresponding author. Mailing address: Department of Therapeu- 
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mM spermidine-HCl (pH 7.5) and were lysed by the addition 
of 300 p.g of lysozyme per ml. The lysed cells were sonicated 
briefly to reduce the viscosity before pelleting the cellular 
debris. The supernatant was loaded directly onto a DEAE 
Biogel column (4 by 30 cm), and the column was eluted with 1 
liter of a linear gradient of salt (0 to 500 mM NaCI). Fractions 
(6 ml) were collected and assayed for TK activity (see below). 
Fractions containing TK activity were pooled and applied to a 
hydroxylapatite column (1 by 8 cm) equilibrated with 1 mM 
potassium phosphate buffer (pH 6,8). The column was washed 
in succession with one column volume of the equilibrating 
buffer, the equilibrating buffer containing 500 mM NaCI, and 
finally the equilibrating buffer. The column was eluted with a 
100-ml linear gradient of 10 to 300 mM potassium phosphate 
(pH 6.8) (10). Fractions with TK activity were dialyzed against 
0.2 mM thymidine-50% glycerols mM Tris (pH 6.8) and 
were stored at -20°C. The enzyme was judged to be about 
30% pure on the basis of the staining intensity of the TK 
polypeptide on polyacrylamide gels. 

Standard TK assay. Enzyme activity was determined by the 
differential binding of phosphorylated versus unphosphory- 
lated nucleoside to positively charged DEAE paper (26). All 
assays were stopped prior to the consumption of greater than 
5% of the substrate to ensure the linearity of the assay. 

Assay for inhibition by competition. To assess the binding of 
nucleoside analogs to the TK, cold nucleoside analogs were 
used to compete with [ 14 C]thymidine as a substrate. To the 
reaction mixture was added 2 mM unlabeled nucleoside ana- 
log. The concentrations of the other reagents were kept the 
same. The assay conditions were identical to those described 
above for the standard TK assay; in particular, thymidine was 
present at 0.2 mM. 

P transfer assay. Phosphorylation of the nucleoside analogs 
was determined by the transfer of radioactive phosphate from 
[•y- 32 3P]ATP to the nucleoside. This assay was performed as 
described by Doberson and Greer (8). The specific activity of 
[ 32 P]ATP used in these assays was 18 |iCi/mmol. 

The open reading frame from nucleotides 144861 to 143038 
on the EBV genome (BXLFI) has previously been identified as 
having some homology to the HSV-1 and HSV-2 TK genes 
and, hence, is likely to encode an EBV TK (3). As predicted, 
the expression of this open reading frame yielded a polypep- 
tide of about 70 kDa. This protein was found to aggregate in 
concentrated solutions and in low-ionic-strength buffers, mak- 
ing purification difficult. In this respect the recombinant en- 
zyme behaved similarly to the TK activity from induced, 
EBV-containing lymphoblastoid cells (7). Gel filtration of the 
soluble form supported the conclusion that the enzyme exists 
in solution as a monomer of the expected size. 
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FIG, 1. MetaJ ion usage by EBV TK. No activity was observed in 
the absence of metal ions. Assay conditions were 10 mU metal 
chloride, 0.2 mM [ 14 C)thymidine (0,05 Ci/mM), 10 mM ATP, 0.6 mg of 
BSA per ml, and 100 mM MES (pH 6.0). 
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Initial characterization of the EBV TK activity was directed 
at determining the divalent metal ion requirements and the 
optimal pH for TK activity. EBV TK can use a variety of 
metals, with magnesium providing optimal activity (Fig. 1). 
Copper and iron are ineffective, while barium, calcium, cobalt, 
and manganese provide roughly the same activities. For sub- 
sequent study of substrate affinity, magnesium was the divalent 
ion used. 

The apparent affinities of the EBV TK for the thymidine and 
ATP substrates were determined by measuring the reaction 
velocity as a function of each substrate concentration in the 
presence of an excess of the other. Lineweaver-Burk plots for 
both thymidine and ATP were linear (Fig. 2), This result 
suggests that the enzyme can be described by the simple 
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FIG. 2. Lineweaver-Burk plots varying dT and ATP. The K m values 
obtained were 22 |iM for dT and 25 u.M for ATP. 



Percent Inhibition of TK Activity 

FIG. 3. Inhibition of EBV TK activity by nucleoside analogs mea- 
sured by competition with thymidine. Assay conditions are described in 
the text. 



Michaelis-Menten model for enzyme behavior. The apparent 
K m for thymidine is 22 u,M, and that for ATP is 25 p,M. 

Competition assays between nucleoside analogs and radio- 
labeled thymidine showed that a variety of nucleoside analogs 
can compete with thymidine for EBV TK (Fig. 3). Several 
5-halogenated deoxyuridine analogs (5-chlorodeoxyuridine 
[CldU], BrdU, and 5-iododeoxyuridine [IdU]) appear to be 
even more effective competitors than the normal substrate 
thymidine. Although these comparisons are not as quantita- 
tively precise as complete KJKi determinations, they allow 
preliminary estimations of the potential usefulness of a given 
analog as an alternate substrate by the EBV TK. Since the 
molar ratio of the test compound to [ 14 C]thymidine was 10:1, 
this experimental approach is most informative for inhibitors 
which compete well with thymidine. 

To test if the competitive analogs were able to function as 
alternate substrates and not just as inhibitors of true substrate 
binding, we examined the formation of phosphbrylated prod- 
uct by using the appearance of hew 32 P-labeled material as the 
basis for such activity. This assay is able to detect agents which 
are relatively efficiently phosphorylated; very low levels of 
phosphorylation may be masked by the unavoidable back- 
ground of radioactive material which appears even in the 
absence of enzymatic activity. The following compounds exhib- 
ited detectable competition at concentrations 10-fold greater 
than the [ 14 C]thymidine concentration: 5'-aminodeoxyuridine 
(aminodU), 5-bromodeoxycytidine (BrdC), deoxyuridine (dU), 
methyluridine (meU), 4-thiodeoxyuridine (thiodU), 5-hydroxy- 
methyldeoxyuridine (HOmedU), 5-iododeoxycytidine (IdC), 
thymine arabinofuranoside (araT), dideoxythymidine (ddT), 
3'-<9-acetyIthymidine, (acetyldT), 5-fluorodeoxyuridine (FdU), 
5-bromovinyldeoxyuridine (BrVdU), zidovudine (AZT), 5-eth- 
yldeoxyuridme (EtdU), CldU, BrdU, and IdU (Fig. 3). All of 
these competitive nucleoside analogs were also phosphory- 
lated by EBV TK The following nucleoside analogs were not 
recognized by EBV TK: deoxyribosyladenine (dA), deoxyribo- 
sylcytidine (dC), deoxyribosylguanine (dG), a-deoxythymidine 
(adT), 2-deoxyribose, 5'-aminodeoxythymidine, azadeoxycyti- 
dine, dihydroxypropylmethylguanine (DHPG), acyclovir, and 



Vol, 38, 1994 



NOTES 2177 



TABLE 1. Comparisons of dA, dC, and dG with thymidine as 
competitors and substrates for the EBV TK" 



cpm in: 

Nucleoside — 

Competition assay Phosphorylation assay 



None 


1,002 


890 


dT 


190 


3,777 


dA 


1,053 
1,259 


636 


dC 


647 


dG 


1,086 


623 



" For the competition assay each nucleoside was added to the standard assay 
reaction mixture at a 10-fold molar excess over the ,4 C-labeled thymidine 
substrate. For the phosphorylation assay, the production of phosphorylated 
nucleoside was assayed by the amount of 33 P transferred from ATP to the 
indicated nucleoside. 



dideoxycytidine. They were neither competitive with thymidine 
nor phosphorylated. The data for dA, dC, and dG are pre- 
sented in Table 1, 

The substrate behavior of BrdU was analyzed in more detail 
(Fig. 4). [ 3 H]BrdU was used in place of [ ^CJthymidine to assay 
enzyme activity, and these results confirm the fact that EBV 
TK phosphorylates BrdU. With a measured K m of 15 p,M for 
BrdU, EBV TK may have a slightly higher affinity for the 
nucleoside analog BrdU than ifor its presumptive natural 
substrate thymidine. 

TK provides the pathway for phosphorylating various nucle- 
oside analogs important in the treatment of herpes simplex 
virus (6) and human immunodeficiency virus. We character- 
ized some aspects of EBV TK with an emphasis on its ability to 
use various nucleoside analogs. The requirement for a divalent 
metal ion is typical of most ATP-requiring enzymes since most 
such enzymes recognize Mg 2+ -ATP as the physiological sub- 
strate (18). Since, except for Mg 24 *, only one concentration of 
metal ion was examined in our survey, we cannot be sure that 
we have determined the optimal activity; we can be sure only 
that the metals tested function to some extent. 

Linear double-reciprocal plots were obtained for both thy- 
midine and ATP. This result indicates that EBV TK exhibits 
Michaelis-Menten kinetic behavior. For EBV TK the K m for 
thymidine is 22 jaM, while for HSV-1 TK the K m is in the 0.5 
fjtM range (5). Conversely, ATP binding is tighter for EBV TK, 
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FIG. 4. Lineweaver-Burk plot varying BrdU. The K m value ob- 
tained was 15 p.M. Assay conditions were 10 mM NaF, 10 mM MgCl 2 , 
10 mM ATP, 0.6 mg of BSA per ml, 100 mM MES (pH 6.0), and 
various concentrations of [ 3 H]BrdU (0.1 to 0.2 Ci/mM). Data from two 
separate experiments are plotted as open and closed circles. 



with a K m of 25 u,M, as opposed to a K m of greater than 100 
U/M for HSV-1 TK. Since the intracellular poo! concentration 
of thymidine is lower (21) than the K m of the EBV TK, this 
pathway probably is of minor importance except during times 
of lytic viral replication when the cellular DNA may be 
degraded, raising the thymidine concentration. 

It is interesting that, on the basis of the amino acid sequence 
homology, EBV TK and HSV TK have more homologous 
ATP-binding sites than thymidine-binding sites (3, 17). HSV 
TK is known to function also as a thymidylate kinase (4), and 
it will be interesting to see in further studies if EBV TK, 
purified to homogeneity to eliminate the possibility of contam- 
inating proteins, has such an activity. 

As shown in Fig. 3 it is apparent that a variety of nucleoside 
analogs are used by EBV TK. Changes in the base or changes 
in the sugar of the nucleoside can be accommodated by this 
TK. On the basis of the competitive activities of these analogs, 
we can infer that there are regions of the nucleosides in which 
structural modifications can be tolerated without destroying 
substrate-binding activity. Substitution on the sugar is toler- 
ated at the 3' position, with the bulkiness and charge of the 
group possibly determining the acceptability by the enzyme 
(compare AZT with acetyldT), Interestingly, ddT with a hy- 
drogen substituted at the 3' position is not as effective a 
substrate as AZT or acetyldT, which have much bulkier 
constituents. Only the substituents azide and acetyl, however, 
can function as hydrogen bond acceptors. This H bonding may 
be more important in TK substrate recognition than steric 
occupation. 

The 2' position of the sugar may also be substituted, 
although with limited acceptance by the TK. The ribose of 
methyl uridine and the arabinose of araT were both phosphor- 
ylated. Limited phosphorylation of araT by EBV TK indicates 
that the enzyme can distinguish the '2' hydroxyl group in the 
beta orientation of the arabinose sugar. 

Large changes in the nucleobase were least well tolerated. 
EBV TK only binds and phosphorylates pyrimidines. Of the 
naturally occurring nucleosides, thymidine is the best sub- 
strate, with some activity seen with deoxyuridine and no 
activity observed for dC. TK is thus dissimilar to HSV TK, 
which can phosphorylate most pyrimidines and some purines. 
However, Stinchcombe and Clough (23) reported that an 
activity in a DEAE column fraction, which seemed to be 
specific to lymphoblastoid cell lines expressing EBV lytic cycle 
genes, phosphorylated dC as efficiently as it did dT. Likewise, 
Littler and colleagues (17) have reported the ability of EBV 
TK to phosphorylate both deoxypyrimidines and deoxypurines. 
The purine analogs acyclovir and DHPG were also phosphor- 
ylated at a very low level in the preparation of EBV TK 
described by Littler et al. (17). Since the enzyme used by 
Stinchcombe and Clough (23) was an. impure preparation 
isolated from lymphoblastoid cells in culture, it may have 
contained several different activities, some of which were 
unrelated to EBV yet which may have modulated the substrate 
behavior of the EBV TK polypeptide or altered its stability in 
some way. Although our observations, too, are based on only 
partially purified EBV TK, such other contaminating activities 
may have been eliminated from our preparation or may be 
absent from material produced from recombinant genes in E. 
coli. Another possibility to explain the discrepancy between 
our results and those of Stinchcombe and Clough (23) may be 
the phosphorylation of dC by a virus-induced kinase analogous 
to the cytomegalovirus protein kinase which can use dC as an 
alternate substrate (16, 24). 

Deoxyuridines substituted at the 5 position were effectively 
used by the enzyme. Substituents at this position ranged from 
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vinyl to halogenated compounds. The chloro, bromo, and iodp 
5-substituted deoxyuridines were the most effective nucleoside 
analogs and were better even than the physiological substrate 
thymidine. A comparison of the van der Waals radii suggests 
that these halogens (CI = 0.18 nm, Br « 0.195 ran, I = 0.215 
nm) are very close in size to the methyl of thymidine (0,20 ran), 
On the basis of steric factors, these analogs should bind as 
effectively as thymidine. With a somewhat smaller size, FdU 
(0.135 nm) was the least effective competitor among this series 
of 5~halodeoxyuridine analogs. However, the bromovinyl sub- 
stituent is tolerated despite its larger size. The terminal 
bromomethene group on BrVdU may be oriented away from 
the active site and the proximal methene group may be situated 
to interact with the active-site pocket. This orientation would 
make BrVdU sterically similar to dT. A similar argument can 
be made for the proximal methene group of EtdU. 

Hydrogen bonding at the active site may also play an 
important role in substrate recognition. Because of the elec- 
tronegativities of the halogens, their effects on the nucleobase 
are to favor the enol tautomer as opposed to the keto tautomer 
favored by all other nucleobases, including thymidine (9). Such 
a change in hydrogen-binding potential, along with steric 
factors, may provide augmented binding in the active site and 
result in a substrate with a higher affinity than that of the 
naturally occurring one. That electronegativity alone is insuf- 
ficient to explain pur results, however, was shown in the case of 
the fluoro derivative, which, although it was the most electro- 
negative, was a poor substrate. The steric effect in the case of 
the fluoropyrimidines may negate the increase in H bonding 
potential. 

ACV has been reported to inhibit EBV replication in cell 
culture (13, 15) and has been used to treat EBV infections in 
a number of different clinical settings (1, 11). Our observations 
of the inability of EBV TK to phosphorylate this nucleoside 
analog are not incompatible with these findings. Although 
several studies have shown that ACV is inefficiently phosphor- 
ylated by cellular or virus-induced enzymes (14, 22), the action 
oi ACV as a DNA chain terminator means that even very low 
levels of phosphorylation can be effective in stopping DNA 
replication. This interpretation is compatible with the litera- 
ture reports of ACV efficiently inhibiting the production of 
EBV in culture and human patients whUe at the same time 
being poorly phosphorylated. While the clinical use of ACV 
has shown some success in the case of oral hairy leukoplakia, 
an EBV-related condition, other sorts of EBV infections have 
been reported to be resistant to treatment with ACV (2, 19, 
27). 

In summary, we characterized some aspects of EBV TK 
function, including the use of various nucleoside analogs to 
functionally probe the enzyme's active site. Understanding of 
the active site of this enzyme will allow a better understanding 
of the enzyme's mechanism of action and, in turn, the design 
of better substrate analogs with possible therapeutic impor- 
tance. 

This work was partially supported by NIH grants POl-GM 39546 
and POl-CA 39238. Peter P. Tung was partially supported by NIH 
postdoctoral training grant T32-09159. 
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Abstract 

Background: Development of countermeasures to bioterrorist threats such as those posed by 
the smallpox virus (variola), include vaccination and drug development. Selective activation of 
nucleoside analogues by virus-encoded thymidine (dThd) kinases (TK) represents one of the most 
successful strategies for antiviral chemotherapy as demonstrated for anti-herpes drugs. Vaccinia 
virus TK is a close orthologue of variola TK but also shares a relatively high sequence identity to 
human type 2 TK (hTK), thus achieving drug selectivity relative to the host enzyme is challenging, 

Results: In order to identify any differences compared to hTK that may be exploitable in drug 
design, we have determined the crystal structure of WTK, in complex with thymidine 5- 
triphosphate (dTTP). Although most of the active site residues are conserved between hTK and 
WTK, we observe a difference in conformation of residues Asp~43 and Arg-45. The equivalent 
residues in hTK hydrogen bond to dTTP, whereas in subunit D of WTK, Asp-43 and Arg-45 adopt 
a different conformation preventing interaction with this nucleotide. Asp-43 and Arg-45 are 
present in a flexible loop, which is disordered in subunits A, B and C. The observed difference in 
conformation and flexibility may also explain the ability of WTK to phosphorylate (South)- 
methanocarbathymine whereas, in contrast, no substrate activity with hTK is reported for this 
compound. 

Conclusion: The difference in conformation for Asp-43 and Arg-45 could thus be used in drug 
design to generate WTK/Variola TK-seiective nucleoside analogue substrates and/or inhibitors 
that, have lower affinity for hTK. 



Background 

Thymidine kinases form part of the salvage pathway for 
pyrimidine deoxyribonucleotide synthesis, TKs are 
expressed in a variety of organisms from human to bacte- 
ria as well as in a number of viruses. The reaction catalysed 



by TK involves the transfer of a y-phosphoryl moiety from 
ATP to 2'deoxy-thymidine (dThd) to produce thymidine 
5'. -monophosphate (dTMP). Certain TKs, such as those 
from herpes simplex virus type 1 (HSV-1) and varicella 
zoster virus (VZV) have, in addition, thymidylate kinase 
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activity allowing the conversion of dTMP to thymidine 5'- 
diphosphate (dTDP). TKs can be classified into two types 
which differ in several respects [1]. Type 1 TKs are of 
higher molecular weight, typically around 40 kDa, and are 
active as hornodimers. This subfamily contains the HSV1, 
HSV2 and VZV TKs, and also mitochondrial TK. 

TKs of type 2 include those from poxviridae such as vac- 
cinia virus (W) and variola virus, [2], as well as from 
human [3] hTK, (human type II thymidine kinase 1) and 
mouse [4], Type 2 TKs have a smaller polypeptide chain 
compared to type 1, being of ~25 KDa but form homote- 
tramers. They are sensitive to the feedback inhibitors 
dTDP or dTTP, which are generated at the end of the met- 
abolic pathway [5). Type 2 TKs have a much narrower sub- 
strate specificity compared to type 1 TKs and only 
phosphorylate 2'deoxyuridine (dU) and/or dThd [6]. For 
example, the antiherpetic drug (E)-5-(2-bromovinyl)- 
dUrd (BVDU) [7] is not metabolised by the type 2TKs (i.e. 
cytosolicTK) in contrast to the type 1 TKs (i.e. mitochon- 
drial TK, HSV-1 TK) which can phosphorylate a variety of 
(5-substituted) nucleoside analogues including BVDU. 
Moreover, HSV-1 and HSV-2 TK can even recognize (acy- 
clic) purine nucleoside analogues such as acyclovir and 
ganciclovir [8]. This difference in substrate specificity is 
the basis of some selective antiviral drug action as vali- 
dated by the activation of nucleoside analogues by certain 
herpes virus TKs. Moreover, herpes TKs are also being 
studied as suicide genes in a combined gene/chemother- 
apy strategy to treat cancer [9]. 

The World Health Organisation declared in 1980 that 
smallpox had been eradicated. Since then, routine inocu- 
lation with the vaccinia virus vaccine was discontinued, 
resulting in minimal or even non-existent smallpox 
immunity in the human population [20]. Today, the 
potential use of smallpox virus as a biological weapon is a 
major cause for concern, particularly in the context of cur- 
rent low levels of herd immunity to the virus. Addition- 
ally, the re-emergence of monkeypox virus infection in 
humans (mainly in Africa but some cases have also been 
reported in the United States [11]), has lead to the stock- 
piling of smallpox vaccine (W), mainly in developed 
countries [12]. Nevertheless, some adverse reactions 
which are sometimes lethal, following vaccination have 
been reported [13-15]. W should neither be given to preg- 
nant women for example, nor to people who have a weak- 
ened immune system, skin problems like eczema, heart 
problems or to children under one year old [12). Thus, 
specific anti-variola drugs need to be developed as a mat- 
ter of priority, particularly for widespread use in a bioter- 
rorism emergency, as well as for specific cases of 
unwanted contamination by W or complications like 
eczema vaccinatum or progressive vaccinia. Such drugs 
would be of particular importance if more virulent strains 



of variola virus were engineered using genetic modifica- 
tions, Worryingly, one of the only anti-variola drugs avail- 
able, cidofovir [16], has been reported to be inefficient 
against several pox viral strains [17]. 

Herpes virus TKs have been exploited in the selective acti- 
vation of nucleoside analogues such as acyclovir, which 
have spawned a series of highly successful anti-herpetic 
drugs [18). The same approach may therefore be applica- 
ble in the case of orthopoxviruses. Further understanding 
of the structural differences between type 1 and 2 TKs 
would also greatly help to improve and/or create specific 
drugs against the type 1 TK. Indeed, a drug would be selec- 
tive if it is preferentially metabolised by the type 1 and not 
the type 2 TKs or vice versa. Thus gathering structural infor- 
mation about both TK types is likely to be of great impor- 
tance in assisting drug design. So far, HSV1 TK [19,20] and 
VZVTK [21] structures have been solved, whilst for type 2 
TK, the structure of human cytosoiic thymidine kinase 
[22,23] and the TK of Ureaplasma ureolyticum recently 
became available [22,24], Moreover, determining the 
structure of different type 2 TKs could help as in the case 
of orthopox viruses, to design virus TK-specific drugs. 
Recently, 5-(2-amino-3-cyano-5-oxo-5,6 / 7,8-tetrahydro- 
4H-chromen-4-yl)-l-(2-deoxy-pentofuranosyl)pyrimi- 
dine-2,4-(lH,3H)-dione has been shown to be phophor- 
ylated specifically by vaccinia and cowpox virus TKs [25]. 

We report here for the first time, the structure of WTK, 
one of the smallest type 2 TKs known, in complex with 
dTTP at 3.1 A resolution. This work will be of use, in com- 
bination with the previous type 2 TK structures, to design 
or improve type 1 TK-selective drugs as well as drugs that 
show selectivity against orthopox virus TKs such as WTK 
or Variola TK. 

Results and discussion 
WTK enzyme kinetics 

Assessment of the activity of WTK to phosphorylate vari- 
ous nucleosides showed its ability to recognize .both dCyd 
and dThd amongst the natural substrates. The efficiency of 
dCyd phosphorylation was, however, much less (<5%) 
compared to that with dThd (Fig. 1). In contrast, purified 
cytosoiic TK showed a poor, if any, recognition of dCyd 
(data not shown). WTK also efficiently converted the thy- 
midine analogues araT and AZT to their corresponding 
mono-phosphates whereas no substrate activity was 
observed for the purine nucleosides araA, dA, dG and 
CdA. In this respect, WTK behaves more similarly to 
cytosolicTK (Fig. 1). Also, dTTP inhibits the phosphoryla- 
tion of dThd (1 uM) by WTK with an IC S0 .(50% inhibi- 
tory concentration) of 14 ± 4.0 uM, whereas the 
corresponding IC 50 for cytosoiic TK was markedly lower 
(2.8 ± 0.5 fiM). 
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Figure I 

Histogram showing substrate activity of certain nucleosides 
with WTK. 



Overall structure of WTK 

The quaternary structure of WTK (Fig. 2a and 2b) is tetra- 
meric which is similar to that of hTK [22,23]. hTK shares 
66% amino acid sequence identity with WTK, showing a 
conserved fold amongst species (Fig. 2b). The four WTK 
subunits within the tetramer are identical apart from a 
flexible loop (residues 40 to 60) situated on the surface. In 
the other type 2 TK structures [22-24], the equivalent loop 
is either not visible or appears to be in a different confor- 
mation from one monomer to the other, indicating that 
this flexibility is common to Type 2 TKs. It has been 
inferred that this loop might be involved in phosphoryl 
donor interaction [24]. 

As found in other type 2 TKs, WTK has a central oc/P 
domain structure consisting of six parallel p-strands sur- 
rounded by helices (Fig, 2a). Above the a/p domain, a 
zinc-binding domain is present, the latter belongs to the 
structural zinc-binding class [26]. Although no zinc was 
added during purification or crystallization, the metal is 
present in the structure and coordinated to four cysteines. 
This metal is believed to play a structural role, stabilising 
the adjacent loop, or lasso domain, involved in nucleo- 
side binding via residues Tyr-166 and Arg-150. The 
cysteines involved in zinc binding are also found in the 
hTK, whilst in the case of the Ureaplasma enzyme, the last 
cysteine is substituted by a histidine [22]. 

Active site structure of WTK 

WTK crystals grew in the presence of dTTP, the feedback 
inhibitor of the pathway. The binding of dTTP in WTK is 
similar to that observed in the hTK structure. The phos- 
phates of dTTP bind to the ct/P domain, whereas the 
thymine base and the deoxyribose bind between the lasso 
domain and the a/p domain (Figs. 2a, 3c). The thymine 



moiety is bound to the backbone of the protein via resi- 
dues Phe-113, Ile-157 and 159. Thus this region of the 
enzyme appears unsuitable for use in designing species- 
selective ligands as no side-chain contacts are present (Fig. 
2c). 

Whilst most of the WTK active site residues are conserved 
and superimpose well with the human TK, nevertheless, 
some significant differences are apparent. In the WTK 
active site Ser-148 replaces the equivalent residue in hTK, 
Thr-163, a serine being also present in UuTK. Welin et ah 
proposed that this position, relatively close to the 5- 
methyl group of dTTP, could be used in the design of 
selective nucleoside analogues for UuTK versus the 
human enzyme. Indeed the pocket around the 5 -methyl 
group of dTTP is slightly larger compared to the pocket of 
hTK, so substitutions at the 5 position of dThd could be 
explored [24]. Based on the same argument, these kinds of 
analogues may also be effective with WTK and UuTK but 
not hTK. 

A further important difference in the WTK active site is 
the conformation of Asp-43 and Arg-45 (Fig. 3b). The 
equivalent residues in hTK are in close contact with the 3'- 
oxygen of the deoxyribose and to the P-phosphoryl oxy- 
gen of dTTP respectively, whereas they are not . in the case 
of WTK. Asp-43 and Arg-45 are only observed in the 
WTK subunit D, they are disordered in the other three 
subunits, indicating the flexible nature of this loop. Nev- 
ertheless, the dTTP still adopts the same conformation in 
each subunit as observed in the human enzyme. We thus 
can infer that Asp-43 and Arg-45 may not be crucial for 
the binding of dTTP, otherwise this ligand would either 
not stay bound to the active site or would shift position. It 
has been proposed that Arg-60 in hTK (Arg-45 in WTK) 
may act to stabilize the transition state. of the reaction 
[22]. From the WTK structure, it appears that the confor- 
mation of this residue has to change, to allow contact with 
the substrate and hence the reaction to occur. The flexibil- 
ity of this loop region also implies that more space is 
available in the WTK active site compared to hTK. Clearly, 
using bulkier groups at the 3- and 5-positions in the pyri- 
midine ring of the dThd molecule should be investigated 
for the rational design of novel, specific drugs against 
WTK. Modifying the sugar moiety would also be a worth- 
while approach to pursue in the design of further ana- 
logues. 

Ligand specificity of WTK 

Thymidine analogues such as (South)-methanocarbath- 
ymine ((S)-MCT) and (North)-methanocarbathymine 
((N)-MCT) (Fig. 4) have been studied as antiviral agents 
against orthopox and herpes viruses. For orthopox viruses, 
such as variola and vaccinia viruses, ((N)-MCT has potent 
antiviral activity whilst ((S)-MCT) does not [27,28], Due 
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Figure 2 

a. Diagram showing one WTK monomer complexed with dTTP and magnesium, a-helices are shown in red/grey, (3-strands 
are in blue/grey with loops in gold. dTTP is in standard atom colors and the magnesium ion is in blue. b. Diagram showing the 
superimposition of a-carbons of hTK (grey) and WTK. Each WTK subunit is shown in different colors (green/red/yellow/ 
blue), c. Stereo-diagram of the WTK active site showing residues that interact with dTTP and magnesium. Protein side-chains 
and main-chain carbon atoms are in grey, main-chain N and O atoms are in blue and red respectively). dTTP is drawn in stand- 
ard atom colors with the magnesium ion in green 



Page 4 of 9 

(page number not for citation purposes) 



BMC Structural Biology 2006, 6:22 



http://www.biomedcentral.eom/1472-6807/6/22 





C 




a. Stereo-diagram highlighting the positions ofresidues studied in site-directed mutagenesis experiments of WTK. Side chains 
subjected to mutation are shown in grey. b. Stereo-diagram showing the different conformations of residues Asp-43 and Arg- 
45 in subunit D of WTK and the equivalent hTK residues colored in cyan. c. Diagram showing the final 2fo~fc map contored at 
la for dTTP and magnesium. Color coding for the protein secondary structure is as for Fig 2a. 
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to the pseudo-ring conformation, ((S)-MCT) has been 
reported to be the preferred substrate of herpes TKs, 
whereas ((N)-MCT) is preferred by the cellular DNA 
polymerases [29]. Prichard et al suggested that ((N)-MCT) 
could also be phosphorylated by type 2 TKs from cowpox 
and vaccinia viruses [28], whereas Smee et al inferred that 
orthopoxvirus TKs are not responsible for the formation 
of phosphorylated ((N)-MCT) [27]. it has also been 
shown that hTK has a weak affinity for ((N)-MCT) and 
essentially no affinity for ((S)-MCT) [29]. ((S)-MCT) 
inhibits growth of HSV1 TK transduced osteosarcoma 
cells [30], thus hTK alone could not be involved in drug 
activation. The MCT activity seems to be highly dependent 
on the particular cell line, which is thus a consideration in 
validating a particular result, Thus the test system ulti- 
mately of most relevance to human therapy would 
involve primates [27]. 

Nevertheless, HPLC activity assays have shown that ((S)~ 
MCT) can be phosphorylated by WTK [23], thus MCT 
modeled in the WTK active site structure could be used as 
a basis for drug design. The superimposition of (S)-MCT 
and dTTP shows that the thymine ring and the sugar/ 
pseudo-sugar match very well. As the (S)-MCT can be 
phosphorylated, this suggests that the conformation of 
the dThd moiety of dTTP has a very similar conformation 
to the dThd substrate in WTK. This has been confirmed 
by the recent structure of UuTK complexed with dThd 
(PDB code: 2B8T) [24]. The fact that WTK is able to phos- 
phorylate (S)-MCT suggests that a larger sugar moiety can 
fit into the active site. Such a difference in substrate specif- 
icity between hTK and WTK might be due to the flexibility 
of the loop containing Arg-45 (Arg-60 in hTK). 



Further suggestion that the plasticity of the WTK active 
site relative to hTK is indicated by the selective phosphor- 
ylatioin by WTK compared to hTK of the bulky molecule 
5-(2-amino-3-cyano-5-oxo-5,6,7,8-tetrahydro-4H- 
chromen-4-yl)-l-(2-deoxy-pentofuranosyl)pyrimidine- 
2,4-(lH,3H)-dione. The latter drug showed a minimal 
toxicity to uninfected cells and inhibited the replication of 
vaccinia and cowpox virus [25]. Fan et ah suggested that 
orthopoxvirus kinases are more promiscuous than previ- 
ously believed. Indeed their substrate specificity seems 
broader than cellular kinases, thus these enzymes could 
be exploited in a similar manner to the herpes TKs. 

Relating WTK mutagenesis and structural studies 

Prior to the determination of any crystal structure of a type 
2 TK, some mutational studies had been reported on 
WTK. It has been shown that residue Asp-82 participates 
in the binding of magnesium [31] (Fig. 3a) and that the 
loss of a negatively charged residue at this position alters 
the ability of WTK to transfer the y-phosphate moiety of 
ATP to dThd. It was suggested that Asp-82 is not involved 
in ATP binding but rather aids its correct orientation for 
nucleophilic attack. In the WTK structure, Asp-82 is 
indeed orientated towards the metal ion and its distance 
from the magnesium ion varies from one subunit to 
another (from 3,2 to 3.7 A). 

WTK mutants Serl8Thr andThrl9Ser (Fig. 3a) [32] have 
also been studied, and showed an increased activity of 
3.7-fold and 1.4-fold, respectively. The equivalent to Thr- 
19 is conserved in most type 2 TKs although not in the 
case of Ureaplasma, where it is substituted by an alanine 
and thus is unlikely to be essential for catalysis. Indeed, 
WTK mutated in this position to the bulkier residues Asn 





Figure 4 

Structure of (north)-methanocarbathyrnine (left) and (south)-methanocarbathvmine (right). 
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or Arg retains significant activity («v27% of wild-type). 
Even a Pro- 19 mutant retains 10% of wild-type WTK 
activity. Thr-19 is not directly involved in dTTP binding 
but may have a role in ATP binding. As yet no type 2 TK 
has been solved with bound ATP, hence it is difficult to 
assess the importance of this residue. The Serl8Thr 
mutant shows a greater enhancement of enzyme activity 
(~3. 8-fold) than does Thrl9Ser and indeed in our struc- 
ture, Ser-18 forms hydrogen bonds via its main-chain 
nitrogen to a -/-phosphate oxygen and via its side-chain 
oxygen to the magnesium and a y-phosphate oxygen. Thus 
Ser-18 would appear more crucial for enzyme catalysis 
than Thr-19. Although Ser-18 is conserved in the hTK 
structure (Ser-33) it is conservatively substituted for thre- 
onine in UuTK (Thr26). The methyl group of the Thr-26 
side-chain is in van der Waals contact with Ile-30 (3.8 A). 
It is thus possible that replacing the serine by threonine 
could stabilize the side-chain oxygen in a preferential con- 
formation for substrate binding, which may explain the 
increased activity of the Serl 8Thr mutant WTK. 

Finally, mutational studies have been reported on Gln- 
1 14 [33] (Fig. 3a). The authors inferred that Gln-1 14 par- 
ticipates in the feedback inhibition by dTTP since this 
mutant TK could bind dTTP yet feedback inhibition was 
abolished. Interestingiy, Gln-114 is not in contact with 
dTTP, and more surprisingly, it is located on the surface of 
the tetramer. The effect of mutating Gln-1 14 to either His 
or Asp would only be to break the hydrogen bond with 
the adj acent Lys . In this case, the abolition of the feedback 
inhibition could only be explained by a destabilization of 
the loop which is in close contact with the dTTP molecule. 
However, the mutation does not apparently affect the 
WTK activity under standard conditions, thus either the 
binding of substrates dThd and ATP are not disrupted by 
this mutation, or if they are, the decrease in apparent 
affinity is compensated for as a result of an increase in cat- 
alytic rate. To clarify this, further experiments would be 
needed to determine the k cal and K m of each mutant. 

Conclusion 

In conclusion, the work reported here provides the first 
structure of WTK. Even though showing similarity to the 
other type 2 TK structures, the detailed comparison of the 
active site points to principles for designing specific inhib- 
itors of WTK relative to hTK, mainly by using the fact that 
the deoxyribose appears able to accept bulkier/modified 
substituents. This result is consistent with Birringer et aVs 
report of (S)-MCT being phosphorylated by WTK whilst 
not being a substrate for hTK [23]. Knowledge of the resi- 
dues involved in ligand binding of type 2 TKs by structure 
comparison, might also assist in the design of specific 
drugs against type 1 TKs. 



Methods 
Protein purification 

The pET-16b-WTK (derived from Novagen pET-lGb) 
expression plasmid coding for an N-terminal His 6 -tagged- 
WTK was transformed into Rosetta(DE3)pLysS for pro- 
tein expression. The cells were grown in a 40 ml Luria 
Broth starter culture supplemented with 50 ug/ml carben- 
icillin and 34 ug/ml chloramphenicol, overnight, with 
shaking at 37 "C. The culture was then diluted 1/100 into 
4 1 of LB and grown at 37 °C until the A 600 reached 0.7. iso- 
propyl-beta-D-thiogalactopyranoside was added to 1 mM 
and the induction was carried out for 3 hours at the same 
temperature. The cells were harvested by centrifugation 
and the pellets were resuspended in buffer A (50 mM 
sodium inorganic phosphate pH 7.8, 300 mM NaCl). The 
cells were disrupted by sonication and the supernatant, 
following clarification by centrifugation, was applied to a 
5 ml nickel chelating column (Hi-trap chelating Column, 
GE Healthcare) pre-equilibrated with buffer A. WTK was 
eluted with buffer A + 500 mM imidazole. Fractions con- 
taining WTK were applied to a Superdex S200 16/60 gel 
filtration column pre-equilibrated with 50 mM Tris-HCl 
pH 7.8 and 200 mM NaCl. The fractions containing WTK 
were pooled and buffer exchanged against 20 mM Tris- 
HCl pH 8 and finally, concentrated to 22 mg/ml. The 
purification protocol produced on average 25 mg of WTK 
per litre of culture. 

Crystallisation and data collection 

Thymidine triphosphate, dithiothreitol and MgCl 2 were 
added to WTK for initial screening in several Hampton, 
Wizard and Emerald kits in a total of 768 conditions. A 
Cartesian Technologies pipetting robot was used to set up 
100 + 1 00 nl sitting drops in Greiner 96-welI plates which 
were placed in a TAP storage vault [19]. After optimisa- 
tion, crystals grew in 2-7% polyethylene glycol 3350, 5 
mM MgS0 4 and 50 mM MES pH6.5, Even though the 
crystals were small (50 um in the longest dimension), 
they were suitable for data collection. 

X-ray diffraction data were collected at 100 K in-house 
using a Rigaku MicroMax 007 generator Cu Ka radiation 
source and a MAR 345 imaging plate detector (MAR 
Research). Crystals diffracted to a resolution limit of 3.1 A. 
Images were indexed and integrated with DENZO, whilst 
data were merged using SCALEPACK [34]. The crystals 
belong to space group P3 E with unit cell lengths of a » b = 
63.3, c = 166.4 (A), and contained four. WTK subunits 
within the asymmetric unit. Detailed statistics for X-ray 
data collection are given in Table 1. 

The WTK crystal structure was solved by molecular 
replacement using MOLREP [35] with the coordinates of 
hTK (PDB code 1W4R) [23] as the search model. Refine- 
ment was carried out using cycles of manual rebuilding 
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Table I: Statistics for crystal I ograp hie structure determination 


Data collection details: 


Data collection site 


In house 


Detector 


MAR345 


Wavelength (A) 


1 .5418 


Resolution range (A) 


30.0-3.1 (3.21-3.10) 


Number of unique reflections 


13636 


Redundancy 


3.8 (2.8) 


Completeness (%) 


94.8 (87.9) 


Average I/g(I) 


8.5 (2.3) 


Rmerge 3 


0.149 (0.494) 


Refinement statistics: 


Resolution range (A) 


30.0-3.1 (3.18-3.1) 


R-factort(R work /R frBe ) b 


25.7/28.9 


Number of protein atoms 


5206 


Average B factor (A 2 ) 


42.3 


Rms bond length deviation (A) 


0.012 


Rms bond angle deviation (°) 


1.42 



Figures in brackets - outer shell data. *R merge - £|l - <I>|/Z<I>; 
bR factor = £|f 0 - FJ/SF,. 



with COOT [36], followed by NCS restrained refinement 
using CNS [37] and REFMAC5 [38] with TLS [39]. Real 
space electron density averaging applied to the four mole- 
cules in the WTK crystal asymmetric unit was carried out 
using the program GAP (D. I. Stuart, J. M. Grimes & J. M, 
Diprose, unpublished). Density for the dTTP ligand was 
fitted as appropriate and refined. The final refinement 
parameters of R work and R^were 25.7 and 28.9% respec- 
tively (Table 1). Structural superpositions were performed 
with SHP [40], and figures were drawn using PYMOL 
(DeLano, W.L. (2002) DeLano Scientific, San Carlos, CA, 
USA). The coordinates of WTK together with structure 
factors have been deposited in the protein data bank (PDB 
code 2J87). 

Enzyme assays 

Thymidine kinase assays were carried out in a 50 jil reac- 
tion mixture containing : 50 mMTrisJHCl pH7.6, 0. 1 mg/ 
ml BSA, 2.5 mM y- 32 [P] ATP, 5 mM MgCl 2 , 5 mM dithio- 
threitol, and nucleoside analogs. In the assays with dThd, 
2'-deoxycytidine, 2'-deoxyadenine and 2'-deoxyguanos- 
ine, 1 mM of nucleoside analog was used, in the other 
assays, 150 julM of nucleoside was used. The samples were 
incubated for 30 min at 37°C. The substrate specificity of 
the purified enzyme was assayed by thin layer chromatog- 
raphy. The autoradiography was scanned with a "Biorad 
Personal Molecular Image FX" and the was interpreted 
with "Quantity one" software also from Biorad. 

Abbreviations 

TK, thymidine kinase; WTK, vaccinia virus thymidine 
kinase; hTK, human type II thymidine kinase 1; UuTK, 



Ureaplasma urealyticum thymidine kinase; dThd, 2 , -deox~ 
ythymidine; dTTP, thymidine S'-triphosphate; dTD(M)P, 
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Understanding the Binding of 5-Substituted 2-Deoxyuridine Substrates to 
Thymidine Kinase of Herpes Simplex Virus Type-1 

I lan.s Do Winter and Piet I lerdowijn* 
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Thymidine' kinase from IISV--1 (HSV-l TK) is a key enzyme in the metabolic activation of 
nniiviral nucleosides. J ligh affinity of such compounds for the enzyme is required for efficient 
phosphorylation. In this study, affinity data from a series of 5-substituted 2'-deoxyuridine 
subsumes in combination with the crystal structure of the viral enzyme were used to investigate 
the structural factors influencing the affinity of these compounds for the enzyme. Calculations 
showed that the binding energetics and conformations of thymidine and the 5-substiiuled 'if- 
uridine analogues are similar. The major part of the binding energy arises from interactions 
involving sugar and base moieties. Small differences in affinity for the enzyme are explained 
by the hydrophobicity of the 5-substituent or by its energetic complementarity with the active 
site pocket. In designing high-affinity nucleoside substrates of HSV-l TK, care should be taken 
lo maintain the geometry of the base moiety and sugar hydroxyl functionalities. Substitutions 
ai the fj position of the nuclcobase should be lipophilic and characterized by well defined 
geometrical properties. The present study represents a first quantitative explanation for 1 ISV-1 
TK affinity of f> substituted 2' deoxyuridincs which are historically the first group of selective 
amivirals" The results may be used to design new and more potent compounds. 



introduc lion 

I lerpes simplex vims {1 ISV) establishes talent infer 
liens in human nerve tissue and exists as type-1 (1 ISV 
!) and type V (I ISV 2). In Tim: lions of HSV-l in humans 
ran lead lo at ule gingivostomat it is, herpes labial is. or 
oculai disease, while 1 ISV 2 accounts for most cases of 
herpes genitalis. 

Antiviral chemotherapy started with the discovery in 
1959 of the aniiherpes activity ol\ r )-iodo-2'-deoxyuridine 
by serendipity. 1 The structure of 5-iodo-2'-deoxy uridine 
is very similar to that of thymidine, i.e. the replacement 
of i he met by! group in t he 5 posit ion of t he uracil moiety 
by a iodine. This .small difference, however, is sufficient, 
lo rentlei a iodo-2'-deo.\yuridine a selective antiviral 
compound. 'I his landmark in antiviral chemotherapy 
was laiei followed by others; the most important are 
( he disco\ en of ai vrlovir as a nontoxic aniiherpes agent 
and the discovers of T a/ido T -deoxyt hymidine as the 
first approved ami HIV nucleoside. The resemblance 
o»" «.h:^ ; a: tei molecule with ! he si ruriurc of i hvmidino 
is also s: i iking. I leu- the Thydroxyl group of the 
deo.wribose mnieiv is replaced b\ an azido group. Hot h 
compounds (a iodo 2' deoxvuridine and [V azido-3' 
deoxyt hymidine) also demonstrate? a vers similar mode 
of action. 1 hey have to he phosphorylaled lo their 
corresponding triphosphate, and in this form they 
selectively intere.it with viral enzymes (UNA poly- 
miM.tse and reverse i ransci iptase. respectively). Ah 
though the stiucime of a-iodo -deoxyuridine and its 
mode ol action have been known for over 30 years, the 
way it interaits with its metabolic and target enzymes 
st ill i emains puzzling. 

Thedisco\ei\ of j iodo deoxyuridine as a selective 
i ISV 1 ill ug lias ied lo the svnthesis of an overwhelming 
number of othei a substituted 2'deoxyuridine com- 
pounds (for a review, see ref X). of which the best known 
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is (/:)-b-(broi no vinyl)- 2' -deoxyuridine/ 1 Structure activ- 
ity studies of this class of compounds have only be cm 
descriptive until now- after the recent discovery of a new 
series of 5-hctoroaromatir-substiiutod 2'dcoxv uridines 
and the availability of the crystal structure of herpes 
simplex type 1 thymidine kinase. 

The genome of HSV-l encodes several virus specific 
enzymes involved in nucleoside metabolism, one of i hem 
being thymidine kinase (TK). This enzyme phosphory 
lates thymidine to its corresponding 5'rnonophosphate. 
using ATP as its phosphoryl source. Therefore. TK is 
a crucial enzyme in the salvage pathway of thymidine- 
5-triphosphate, which is the precursor of the thymidine 
incorporated in DNA. Viral TKs differ in many aspects 
from human TK and are therefore important targets for 
chemotherapy. The viral enzymes have much broader 
substrate specificity than the corresponding human 
enzyme, which only recognizes thymidine as a substrate 
while the viral enzymes even accept [. thymidine, gua- 
nine nucleosides, carbocyclic nucleosides, and many 
5 substituted 2'-dooxy uridine analogues/ 1 

I ISV 1 TK is not essential for viral replication in cell 
culture, but it is an important enzyme for the pathogo 
nicity of the virus. :,h However, also significant foi the 
design of antiviral compounds is that I ISV 1 TK is a 
key enzyme in the metabolic activation of modified 
nucleosides as antiviral agents. The general mode of 
action of these types of antiviral agcrus is that they are 
phosphorylaled to their ^'-monophosphates by ! ISV 1 
TK. The modified nucleosides have a much higher 
affinity for HSV-l TK than for the cellular enzyme, and 
therefore preferential phosphorylation occurs mainly in 
virus-infected cells. The 5' monophosphates are then 
further phosphorylaled to their ^-triphosphates by 
several other cellular kinases, and those triphosphates 
can either act as competitive inhibitors to TIP of the 
viral UNA polymerase or function as substrates for this 
enzyme and become incorporated into the viral DNA. 
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Figure 1. Aiiouiliiu* ol'l ISV 1 thymidine kinase, Thymidine 
CI) and Cln 12!> emphasize the- location of the active site 
porker Loops beivveen residues 149 152. 70-78. and 261 
281 iiif i Missis ig bom I hi 1 crystal structure. Che C- and 
\ urmini .in- labeled (' and N i especlively. 

Thymidine kinase, whose crystal structure has re- 
cently been sulveu at 2.8 and 2.2 A resolution for the 
complexes with ihvinidine and ganciclovir, respectively, 7 
is the first enzyme involved in the metabolic activation 
of antiviial nucleosides. Therefore, it might, be advan- 
tageous lot un effect i\e ami 1 ISV 1 nucleoside to show 
high afliniu Coi I he viral ihvinidine kinase. Thus, any 
project which unlets ihe design of modified nucleosides 
as ant M ISV compounds should incorporate some csti- 
maies on the rale of phosphorylation of the designed 
nucleosides by viral TK. and hence should develop an 
understanding about the interaction of I ISV- 1 TK with 
its subst! aus at the molecular level. 

In ihe present work, t nmputer simulations were used 
u> investigate those interactions and to examine struc- 
tural features important for a subst iluted 2'-deoxyuri- 
dine substrate recognition by I ISV 1 TK. This series 
of a substituted nucleosides, recently discovered by our 
group, an- selectively phosphoryiatod by 1 lSV-1 TK. 8 tl) 
I 'heir admits- foi the enzyme, expressed in terms of lC t( > 
values measured against the natural substrate, ranges 
from to l:iD«M. depending on the substitution pattern 
at the ;i position, these compounds differ from the 
natural substrate thymidine only by a different substi 
union ai the :> position and provide, therefore, an 
interesting, tool for the analysis of the nature of molec- 
ular ivcngniiion between I ISV I I'K and 5-subst ituted 
'A' deoxvuridine analogues. It is expected that the 
results ol this siudv may be of potential interest to drug 
designers and chemists in their search for improved 
i ISV 1 ar.tivirals. Moreover, this study gives a theoreti- 
cal basis lot an historical discovery in antiviral rhemo- 
I hei ap\ . 

The Players 

Ihe Crystal Structure of HSV-1 TK. 1 ISV 1 TK 

has been crystallized in complex with either thymidine 
ni ganciclov ir. 7 The molecule consists of an ab structure 
made up of 13 ix- helices, two 3.,,-holices. and seven 
p'shoets (Figure I}. I he five stranded /i-shoet forms 
pan of the cure ol the molecule, which contains the 
active site Although the sugar equivalents of both 



ligands exhibit similar binding arrangements in the 
active-site pocket, comparison of the nucleobase of 
ganciclovir with that, of thymidine shows that the 
guanine base is flipped by approximately 180°, so that 
its 06 atom is located in a pocket on the opposite side 
of the 5- methyl group of thymidine (Figure 2). As a 
result of these different binding modes, and because of 
the high structural resemblance between thymidine and 
the substrates used for this study, we have taken the 
crystal structure of the thymidine TK complex as the 
starting point for all modeling experiments. It should 
be noted that, apart from the amide function ol' Glu 
1 25 which is rotated by 180°. no other remarkable 
differences between the two structures in the orionta 
lions of the residues lining the pocket are observed 
(Figure 2). 

The crystal structure of HSV-1 TK does not contain 
an ATP molecule. Instead, a sulfate ion is observed at 
the putative position of the /^-phosphate of ATP. 7 This 
sulfate-binding site is located in close proximity to the 
5' hydroxy! function of thymidine (S--05' distance is 6.!i 
A) and is coordinated by the side chains of Arg-222. Arg- 
163, Gln-221. and by the main-chain Nil functions of 
Gly-59, Mel.-60. Cly-61. I.ys-62, and Thr-63. 7 

The Substrates. The 5-substituted 2/-deoxyuridino 
substrates which have been used for the computational 
analysis, together with their affinities for the I ISV- 1 I'K 
enzyme described in the literature, 8 10 are listed in 
Figure 3. Basically, the structures can be divided into 
two classes. 

T he first group involves 5- (thien-2 y!) 2'-deoxvui idine 
(1) as lead compound, with an \O l{i value around 2 «M. 
Affinity for the enzyme is retained when the heterocyclic 
ring is connected via its 3-position (8), when the ring is 
replaced by an isothiazole (7). or when a halogen or 
methyl group is introduced at. the 5 position of the 
r>-(ihien-2-yl) ring {2. 5. and 6), In contrast, halogen 
substitutions at the 3-position of the thion-2-yl ring and 
t he 4 posit ion of t he t.hien-3-yl ring decr ease t he affinity 
by approximately 1 order of magnitude (3. 4. and 9). 

The second group includes 5-(luran-2-yl)-2' dooxvuri - 
dine (10) as lead structure. As in the case of the 
thiophene compounds, affinity is retained when the 
Curan ring is connected via its 3-position (13). I lowever. 
in contrast to the thiophene derivatives, affinity de- 
creases significantly when a halogen group is substi 
luted at the 5-position of the furan-2 yl ring (11 and 
12) or when the heterocyclic ring is replaced by an 
isoxazole (14 and 15). 

Clearly, there is a difference between both classes of 
molecules. Crystal structure analysis revealed that all 
compounds have a quasi-coplanar geometry between the 
fivc-membcred heterocyclic ring and the uracil base 
(Figure 4). 11 M Intramolecular interactions stabilize 
this coplanarity. Interestingly, in the 2-thicnyl series 
of compounds, coplanarity is stabilized by an interaction 
between the OA carbonyl group of the uracil base and 
the sulfur atom of the thiophene ring (Figure 4, struc- 
ture A). A similar orientation for the isot hiaxole dei iva 
live (Figure 4. structure B) has been deduced from its 
NMR spectrum 1 -' and quantum-chemical calculations, 
which predict a stability of 17 k.]/'mol of the N.MR 
observed orientation relative to the other possibility. : 
In contrast, the orientations of the 2-furanyi and is<jx 
azolyl rings are clearly different from the thiophene 
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Figure; 2. SiLMvuvii'w comparing the binding modes of thymidine {dark gray bonds) and ganciclovir (light gray bonds) in the 
active sue pocket of 1 ISV- 1 I K. Oxygen and nitrogen atoms are drawn as light and dark gray spheres, respectively. I lydrogon 
hoi ui.s air mdiriiii'U as clashed lines. There is a remarkable overlap between the amino acid side chains in the two active site 
porkc.s. I'he onlv exception is (;in-l25. whose amide functionality is Hipped approximately 180" between the two crystal forms 
in onlfi in opiimi/.i' hydiogen bonding oppoil unities. 



analogues, as they arc stabilised by intramolecular 
interactions loci wren the oxygen of the note recycles and 
i he 1 l(> atom of t he ui acil base (Figures 4, structures C 
a nil D). : ho oxygons of the furanyl and isoxa/olyl rings 
point tlioi el ore in an opposite direct ion compared lo the 
sulfur atoms oi' the thienyl or isot hiaxolyl rings. The 
s: run tires, of I Ik 1 3 ihionyl and 3-iurfuryl compounds are 
stabilized hv an int i amolerular contact between 04 of 
the uracil base and i 12 of" the hotoroaromatic rings 
( Figure 4. structure K), which is analogous with the 
stabilizing OT-1 hi interaction observed in the crystal 
structure of (/:) :> (bromovinyl)-X' deoxyuridine (h'igure 
*l. structure F). Mi C-are was taken lo maintain these 
specific conformations during the docking and energy 
minimization calculations by the use of appropriate 
t ursiot i angle potent ials. 

the high structural similarity between the 17 5- sub- 
stituted 2' deoxyuridine analogues suggests thai the loss 
of entropy upon binding the enzyme may be similar for 
all compounds, i'he ;i-subst il uted uracil bases are 
planar, conjugated, and rigid systems, making them 
idea IK suited foi structure- affinity studies such as 
i hese 

Results 

Binding of Thymidine. (A) Act ive Site Interac- 
tions. A stereoplot of the crystal structure with thy- 
midine in the act ive site of 1 ISV- ) I K is shown in h'igure 
2.*' i'hroe hydrogen bonds, involving the 3'- and 5'- 
hvdroxyl groups of thymidine, link tin* sugar* moiety of 
thymidine with the enzvme. while the nurleobaso forms 
i wo h\dit>^ic-ii ht)ii(ts thmugh alums \3 and in 
Lie id it ion to t hese hydrogen bond interact ions, t he uracil 
base is sandwiched oi stacked between the < -methyl 
group oi Viet 128 on one side and the aromatic side 
chair, oi I vt 172 along the other side. The 5-methyl 
gniup is located in a small hvdropliobic pocket lined by 



residues of Tyr-132, Ala-167. Ala-168 (not shown), and 
the aliphatic side chain of Arg-163. 

(B) Binding Energetics, A decomposition of the 
total binding energy of the minimized thymidine 
HSV-1 TK complex on a per-residue basis or separated 
into individual contributions for the base and sugar part 
is listed in Table 1. Although the results are slightly 
dependent on the applied dielectric for calculating the 
interaction energies, the main conclusions remain the 
same. More than 80% of the total interaction energy 
comes from the interaction with only six residues, 
namely the four hydrogen bond partners Glu 83. (>)u 
225. Gln-125, Tyr-101, and residues Tyr !72 artd Met 
128. which sandwich the nucleoside base by means of 
slacking interactions. Thus, the major binding forces 
between thymidine and enzyme arise from hydrogen- 
bonding interact ions involving the two hydroxyl groups 
on the sugar moiety and the amide function of the 
nucleobase, together with stacking interactions of the 
base plane. These conclusions are consistent with the 
decreased affinity of either ganciclovir (1C; {J =■ 50/iM). 1 ' 
where the purine base is less well stacked compared to 
thymidine, or acyclovir (1C : ; U > 500 /<M). 17 where only 
one hydroxy] function is available for hydrogen-bonding 
interactions. 

Binding of 5-Subslituted 2' Deoxyuridine Sub 
strates. (A) Exploration of Alternative Binding 
Modes for These Compounds, Considering the struc- 
tural resemblance between iho 5 substituted 2' cleox\ 
uridine analogues and thymidine, t\ binding eonforma 
lion similar to the binding mode of thymidine seems 
reasonable at firsi sight. However, in principle, eight 
different binding mode's are possible on paper, and these 
are shown in h'igure 5. 

The eight possible binding complexes are distin- 
guished either- by the side-chain conformation of Chi 
125 (Table 2; compare binding modes la. Is. 2a, 2s with 
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Figure 3. Tl it- compounds wiih their numbering label and 
oxpei haemal afTiuit\ for I !S\ 1 thymidine kinase 4 . Affinity 
data wore ublained from literature* !i: and inr ihe Lonccntia- 
uuns m|uirrri id inhibit thvmidine phosphoi ylation by 50% 
a»»aiiiM radiolabeled ihvmidine used al a concentration ol' 1 
ff.V). All affinities wvie measured under the same conditions, 
and hi^hh siniilai /\V values {() ?.U u\\) wvn- found in all 



3a 3s la Is), bv the- ot ieni a: ion of the base 1 wiih 
ivsjK'U K; sugar { sw/.^'i//// conformation; la. 2a. 3a. 4a 
ic/sf/.s Is 2s. 3s. 4s). ami by (lit* hydrogen bonding 
pattern between ihe subMrato base and the amide 
fin trt ion u!' C iln IX;) (la. Is. 3a. 3s vs 2a. 2s. 4a. 4s). 

Ilowewr. although in principle 1 eight different ways 
air feasible in dock a subst ifutrd substrates in the 
aeiive site ol' I ISV 1 I K. ii is obvious ibai only a 
small fraction of those il noi only one. will be oner- 
gel u ally feasible ami will represent the "l rue' binding 
mode. I -iin hermore. since 4 it has been shown that the 
f/ hydroxvl fund ions of these substrates are phospho- 
rylated to the same extent and in a similar way as 
the corresponding hydroxyl fund ion of thymidine. 18 t he 
true binding conformation will likely have its 5'-hy- 
droxyl fund ion in a similar oriental ion as found in the 
i t s si at siruetuic of thi' thymidine (and ganciclovir) 
t omplr.V 

hi aide*: to ehtiitla.te ihe most likely binding confor- 
mation horn all the eight possibilities, we have docked 
and eniHgv minimized model compound 5-(t hion -2-y))- 
?.' deow uridine (1 in the active 1 site of 1 ISV 1 TK and 
have subsequent 1 v o\ atuated eaeh of the enght possible 1 
binding me;de\s hase'el on four criteria: (1) the number 
of hvdrogon bonds between subsirate and enzyme. (2) 
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Figure 4. The? crystal structure conformations of several 
S-substiiuled 2' eleoxyui idinc derivatives. 

the root mean squared (ruts) deviat ion of all active site 
residues from the corresponding minimized TK-thymi 
dine positions, (3) the distance between sugar 05' and 
sulfate S atoms, and (4) the interaction energy bet ween 
substrate and enzyme. The results are summarized in 
Table 3. Minimized binding conformations la, 2a, 3a, 
and 4a of substrate 1 in the active site of I ISV 1 TK 
are shown in Figure (>. 

hollowing 1 able 3, we believe thai the most likely or 
true binding conformation resembles binding mode la 
for the following four reasons. First . in bineling confer 
mat ion la we find that ihe distance between the 1 sugar 
QIV and sulfate S is closer than the' corresponding 
distance to t)3'. In fact, this is observed in all eases 
where an unli- orientation between base 1 and sugar is 
adopted (Table 3. binding modes la. 2a. 3a. 4a). if ihe 
sulfate ion is indeed located near the putative binding 
site of the /^-phosphate of ATP. 7 we might expect from 
these? geometries ttiat it is ateau Of)' instead of atom 03' 
which will get phosphorylated by the enzyme. This is 
in agreement with experiment. 18 On the contrary, 
exactly the oppe)site^ is found for the .sw-oriented binding 
conformations (Is, 2s. 3s. 4s), where phosphorylation 
of the 03' atom is expected and therefore deviates from 
experiment tK (Table 3). 

Second, binding mode la allows the? formation of the 
maximum number of hydrogen bonds between substrate 
and enzyme, leading to a much tighter substrate binding 
compared to the other possibilities (Figure 6). For 
example, due to the different orientation of the substrate 1 
in binding modes 2a and 3a compared to la and 4a. 
hydrogen-bonding opportunities of the sugar OS' atom 
with Olu-83 and the sugar 03' with Clu 22b and Tyr 
101 are lost in 2a and 3a. In conformation 2a this loss 
is partially compensated by the? formation of a now 
hydrogen bond between 03' and Clu-83 (Figure G). 
Additionally, while the docked conformation of 4a prior 
to minimization was such that it a!!e)wed for the 
formation of two hydrogen bonds between the 1 substrate 
base and Gin- 125 (Figure 5). minimization however 
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1 All residues wiihui 8 A nC I lie suhsi Kile and I he substrate itself 
uric allowed to lelux. I lie remaining residues wmc kepi con- 
strained •(> iheii rrvsml snm tui'i' ronrdiuaios using a harmonic 
{(nn- ronsianl o!' HHH! k.J'tnol-A ' Nionbondcd interactions were 
udrulaled usinj', a !0 A ruloff distance Kither a distance 
dependent ' /) ni a rnnstnui dii'leci ric (< 1 ) was used far t he 
i al< teal ii» i ai elrct l osial n inlei an ions. "The relat ive com i ibiil ion 
in ihe loial inleraeiion enen»v Ioi each residue is the interaction 
energy hei ween ihe pailicular jesidue and l he whole substrate, 
divided h\ i hi- una! inliT.u'iinn «-ner«;v heiweeu all residues and 
whole suhsi i aie I'hr i elal i\ e cnnirihui ion lo I he lolal interact ion 
euert'V lt)i each suhsi i ale moiel v js l he inlet action energy between 
the paMituini sahsliale muieh and the whole en/vmo, divided 
\i\ [in- Milai inlr: ai imn cneif'v beiwecn I he whole suhsi cai e and 
lin w iinle i :i/\ me. 

disrupts ihr hydrogen bond between the NcZ atom of 
Ciln-lX;") rind ihe base 02 atom, as (here is a tendency 
fur this residue to Hi}) back to its position as found in 
the crystal structure with thymidine (Figure 6). 

Third, taking the natural substrate thymine as a 
reference the leasi distortion is introduced in ihe 
residues lining ihe pocket when substrate is docked in 
a conformation resembling that of la. This is evident 
from the nits deviation between the minimized TK 1 
and TK thymidine complexes, calculated for all resi- 
dues lining the act i\ e si'.e pocket (Table 3). In contrast . 
the binding o! :> substituted X'deoxyuridine substrates 
in ii runfermal iun resembling binding modes 2a or 3a 
hit : oduces se\ ei al slenr clashes with, for example?. Tyr 
!U1 (binding mode 2a) and Arg i G3 (binding modes 2a 
and 3a) (! igure (>). 

I'ourt h. the inlei act ion energy between substrate and 
en/.yme is significanilv more favorable when the sub 
si i ci t t x is docked in binding conformation la compared 
to theoihoi possibilities. This is likely the consequence 
of a more optimal complementarity (both steric and 
electrostatic) between en/.vme and substrate when 
docked in this conformation. In the case of thymidine, 
more* than 80% of the total interaction energy arises due 
to interactions of the sugar hydroxy I groups hydrogen 
bonded to Clu-83. Gin ?2i). and Tyr-JOl. and from the 
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With Gln-125 in a conformation similar as in the thymidine complex:- 



(la) 



(2a) 




Vv 
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(1s) 



With G!n-125 in a conformation similar as in the ganciclovir complex:- 



(3a) 



S 



.o t 



(4a) 



(45) 



figure 5. Fight possible binding orientations to dock a 
o-substituted ^'-deoxyuridine substrate in Ihe active site of 
I1SV-1 TK. All orientations are drawn relative; to same 
fixed reference frame, i c. the orientat ion of the active site 
pocket is assumed to be similar in all eight cases, for a 
definition of -R, see Figure 3. 

Table 2. Classification of the Possible liinding Conformations 
of the [i Suhstituied 1' Deoxyuridine Substrains hy ihe Adopied 
Oriontaiiou of llie Cln \2\) Side Chain-' 

binding mode 



la. Is. 2a. Us 



3.s. 'la. -1s 



C Ca Cji i'y H)2 
Ca C/i Cy C<) . r )B 
Cji Cy C.'rt Xt2 1 1 

(resembling ihe 
values in Ihe 
thymidine njinplcx) 



IVli 
VM) 

(resem!>liiig !!>e 
values in ! he 
l«ai u it lovir complex) 

•■'Listed are ihe approximate slat line, values I'm ihe lorsioi: 
anejes of (mi 1X5 prior lo energy minimi/anon caleulaiions lvalues 
in degrees). 

pyrimidine basc^ which is stacked between Met- 3 28 and 
Tyr-172, and hydrogen bonded to Gin- 125 (Table 1). It 
is therefore no surprise thai an alternative binding 
arrangement lacking one of these interactions may 
result in a significantly decreased interaction energy, 
and hence, affinity. 

(B) Binding Energetics of a Typical 5-SubstiLut 
ed Nucleoside: 5-(Thien-2-yl)-2 -deoxyuridine (1). 
A decomposition of the total binding energy for 5-(ihien- 
2-yl)-2'-deoxyuridine (1) in binding conformation la is 
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Desn ibod in terms of the number of hydrogen bonds between 
substrate and en/vme. i lie rms deviation of all art.ive-sire residues 
from i heir rni responding minimized I K thymidine, positions, the? 
distance between sugm Do' and 1)3' munis and sulfate S atoms, 
and the inte:adion energy bach binding cuiifortiiat ion was 
generated b\ manually docking substrate 1 in the required 
orientation followed b\ enerp minimization until (lie gradient 
dropped below 0.0 1 k 1/mubA. A nonbonded rut oil of 10 A in 
coU|unri ion uil h a dielei n it i onsiaul of t 1 oi ( / was used. 
Ml icsiduos vvuhiu 8 \ of i lit* subsume and the substrate itself 
weie allowed in lela.s \ he n-maining i esidues weie kept ran 
shamed lo :!:e-i n \Mal s! run ore posit ions using a harmonic force 
, tn is.! , it t: ;:! I 000 k J mob \ I ) >e i cpot i : -d niimhei of hydrogen 
nomk 0 111! t nr. dr\ M'.inii !i ms) ,n id distant es between subsirale 
l:\-tiinx\K .mil snlinie S ; < /(( >:>' '• ■ -S) anil f/(()3'-«-S)| an- from I he 
•rmiimi/i-d unmilrws ..iliuialee. vulh <i disl am e dependent di 
r v;i!!t . ' SinaS.a iesuh> v\eie obtained using a constat!! 
dielet i ! :r < I but ai e no! shown lor i lat il v ' The orient at ion of 
the subsuna- base moieu \\\\\\ i e.spe< I to its sugar moiety. In 
the s's // o: ientai ion. ihe base owgeu is positioned above the 
sugm ring, while K)V. :s oitenied awav hum the sugar ring in the 
a////t unloi mat ion 1 \; mil km of hull ugen bonds bet ween subsl rate 
,,-id en/Min 1 ii! ihe minimized complex The definition of a 
hsdrogen bond was so!el\ based on (list an< e criteria, i.e. ihe actual 
tlist a nee bei ween i lie non b\ dr ogen atoms had to be smaller than 
the Mini "I the win del Waals radii id 1.1 A: N. l.li.'i A). J Tltc 
tool mean stjuaied disiame in A between the minimized I K 1 
and ; K ihvmidiue lomplexes calrulaied loi all residues within 
'A A of the se.bsiiale 1 he ihstanie in A between the sugar Of/ 
and sulfa 1 e *s .Hum ' I he distance in A between the sugar ()3' and 
suliau' S .limn flu- miei action eueigy in k.Jmtol biMween 
sabstiate and en/vme ustng i\vn iliffeieni di(rl(M:i.ric cousianl.s {< 
I and ( /) 

o:\imi In : nbh' 1. \'»m \ siniihtr rosulls wore calculated 
ft,] ,dl oiltci tmispoiiml.s ol I i^ir.c 3 (daia not shown). 
As in i he l iisc uf 1 1 ivmidinu. inure than B0'>b ol" tlio tot.al 
binding uiku^n can \n- aittibuiod 10 ihe sugai and 
ntirk'obasc moieties. Strong hydrogon bonds boiwcon 
ihe sugar bydroxvbs and residues Clu-83. Clu-z!X5. and 
J yr - J 0 1 and beiween the amide atoms N3 and 04 of 
the nucleobase and CAr\-\2i). as well as base-stacking 
interactions, are responsible for this tight binding. 

h is cleat thai substilui ion of the 5-inelhyl group of 
thymidine bv a bulkier, unsalurated group such as a 
thiupheue in 1 has no significant effect on the major 
binding interactions involving base and sugar moieties. 
Dorking oi b substituted ?! cieoxyuridino substrates in 
a roulbrmat ion which resembles binding mode la posi 
( ions t he . -suhsi in tied L>roup of the subsl rat es in a small 
pocket uiiure the meih\l group of thymidine is also 
hound. ! he majoi Iraeiion of I ho solvent -accessible 
surface o! i his pocket is ioi med b\ I yi 13^. A moleculai 

j>i aphic s rxiMvise shows t hat nnly a small displacement 
of l yi 13X is needed to erdargc 1 this pocket such that it 
becomes large enough to occupy ^-substituted groups 
like 4 lhien\! moietie.s (bigure 7). 

Understanding the Differences in Affinity of 
5 Substituted 2 -Deoxy uridine Substrates: The 



Role of the 5-Substituted Side Chain. In an attempt 
to rationalize Ihe observed affinity differences between 
the different 5-subsi Ruled substrates, a binding energy 
affinity relationship study was undertaken. The 18 
substrates listed in Figure 3 were alternatively docked 
in the active site of 1 ISV-1 TK y in a binding confornta 
lion similar lo la (Figures 5 and ()). and their orienta 
lions were optimized by means of energy minimization 
techniques. Several energy calculations v\ere performed 
for each substrate, each of them differing in t he applied 
dielectric constant for calculating electrostatic interac- 
tions or in the number of enzyme residues that were 
allowed to relax. Interaction energies between sub- 
strate and enzyme, under the form of nonbonded van 
der Waals and electrostatic terms, were calculated for 
each complex and for each experimental setup to provide 
a measure of the tightness of binding. The results are 
listed in Table 4. 

In order to account for possible differences in hydro- 
phobicily beiween the compounds, free energies of 
solvation for each of the 18 subsl tales were also 
calculated. The electrostatic contribution to the total 
solvation free energy was calculated by solving the 
I'oisson Bolt/mann equation. ,! ' The nonpolar cumpu 
nent of the total solvation free energy was estimated 
front the solvent-accessible surface of each molecule/" 
It was assumed that the differences in solvation free 
energy beiween the 18 compounds could be largely 
attributed to the different substitution patterns at 
position 5. Therefore, to make the computations corn 
putationaliy feasible and to eliminate a possible depen- 
dency of the calculated solvation energy on the confor- 
mation of the sugar moiety, the whole ^'-deoxyuridine 
moiety was replaced by a methyl group. The results 
are listed in Table 5. 

The correlation coefficients for the least squares fits 
between the logarithm of the lQ,r/s and the calculated 
interaction energies, corrected for differences in hydro 
phobicily, are listed in Table 6. Ihe best correlation is 
found for the case where the dielectric constant is 
treated as a distance-dependent function {( - /) and 
where only the residues positioned within 3 A of the 
substrate were* allowed to relax. The least-squares 
equation for this particular case is given by 

ln(lC 50 J ~ 0.1771b; - KCJRT I 51.9. (1) 

in which IF stands for interaction energy in k.J/mol. 
AC'sui is the solvation free energy of the substrate in kj/ 
mol. RT is expressed in k.l/mol. and 1C; 1( , is expressed 
in,«M. Figure 8 shows a plot of this equation. It is clear 
from eq 1 that the affinity for this series of 5-subsiiiuted 
^deoxyuridine substrates can be improved either by 
increasing the side chain's interaction energy (IF more 
negative) or by increasing the side chain's hydrophobic 
ity (AC7 st ,! more positive). 

'The interaction energies in eq 1 have been calculated 
without explicit solvent, using a distance-dependent 
dielectric and by allowing only part of the enzyme to 
minimize. It is therefore important to stress that both 
the slope and intercept of eq 1 have no real physical 
meaning, although in principle some contributions from, 
for example, enzyme desolvation. conformational en- 
tropy, and steric stress should be present. In principle, 
the equation could be expanded to include these factors 
as well. However, we have chosen not to do so. mainly 
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! igure l> > (1 him ^-vl} ?' dcoxvuridine (1) energy minimized in the active site ofHSV-1 'I K. hour of ihe eight possible binding 
conformations are shown (la. 2a, 3a. 4a: Figure 5). Several residues lining the active-site pocket are also shown, ihe energy 
minimized complex beiwern I K and thymidine (dr awn using light -colored bonds) is shown as a reference. 1 he substrate and its 
c«ji i cspuiuling residues are di awn using dark-colored bonds. Hydrogen bonds are displayed as dashed lines. 



because oi dil'firuli ies in accurately estimating them and 
al.su because of the limited number of data lu fit against. 

Disrussiun 

The Binding Conformat ion of 5-Substituted 2'- 
DooxyuHdinc Substrates. On the basis of an explo- 
ration o! ihe ronlormatiunal spare of !)-(lhion-2-yI) T 
deoxvui ithne (1) in the active site ni l ISV 1 TK. we have 
iuuniiilwn ihemosi likeh binding confnnnal inn of those 
kind ni subs' ra!<*s :ost Iv ivsembles dial of \ hvmidine 



(binding conformation la, Figures 5 and G). This is in 
accord with the observation that ihe affinities for ihe 
enzyme of many of those ii-subst iitited deoxyuridino 
substrates (Injure 3) are comparable with die rnrre 
spending value for thymidine, in this binding runlbrma 
lion the sugar hydroxy! groups cue hulrugni bunded in 
C;iu-83. {'Au '/.?.:). and 1 yr 101 . and with ihe pvrimidino 
base slacked between Met-128 and 1 yr \72, and hydro 
gen bonded lo Gln-l^f). Together, these interactions 
represent more than 80% of the total interaction energy 
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Figure 7. t Jnl> :i minor small position shift of Tyr-132 is 
needed U) enlarge the ptwk.fl of the 5-rnethyl group of thymi- 
dine so thai h becomes large enough Lo occupy furano oc 
ihiophene groups. 1 he \ ellow-colored part is the bottom part 
ul the sol vein accessible surface of Die pocket as found in the 
uysial st rue i lire, colored in blue is the new pocket after a small 
position shift of 1 vi 1 32. 
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■All irvihns iiiixau* .it li-.isi one atom positioned within the 
rami! i!:s!.nin at inr sub.Miati wciv allowed m relax. All uihei 
i !•.:«-% writ- Kip: ■ a::M! .mini n> Kieii < i vstal mi urtute roordi 
. i: ,, s :;ssn .. .i -j.i: Inn e . iinsiani ill 1 1)1)1) k I' tool* A '. 



(lanlo 1). and it is thni dure anticipated that an 
alternative binding arrangement lacking one of those 
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" AG',,,,, A6;.| r . and Aoy stands lor nonpolar, liiectric. and total 
solvation free energy, respectively. Values are in kJ/mol. "Cal 
emulated from the solvent accessible surface ai ea (SAS) usin^ A^...:, 
- 4.7 -I 0.040-(SAS). SAS in A* and AC',,,, in kj/rnol. 



Table B. Corrrlaiion Coefficients r for the Least Squares l-'its 
between the l.ogajithui of the ICa/s and the Calculated 
Interaction Knorgies. with the Inclusion of Corri-n ions foi 
Diffeiences in Myclrophobicity beiweeti I he Compounds 

(JA Ta" 3 A 1 A OA 8 A 

"7 T'~ 0 3 i7x " 0 7 t;3 o.: o l 

, ; 0.4 0 2 0.8 0.0 0.3 0.2 



interactions may result in a significantly decreased 
affinity, as is the case fur, for example, acyclovir. 1 ' 

The Role of the 5-Substiluted Heterocyclic Group 
in Fine-Tuning the Affinity. The 5 substituted ?/ 
deoxyuridine analogues differ from the natural sub- 
strate thymidine only by the different substitution 
pattern at the 5-position. The spread in affinity of these 
compounds is narrow and ranges from 0.3 to 143 /vM 
(Figure 3). It is likely that these small differences in 
affinity are the result, of structural differences of the 
5-substiiuted groups, and not from different binding 
conformations of the base or sugar moiety. Any si rue 
fural changes with respect to these moieties would 
result in a much more drastic change in affinity, such 
as in the case of acyclovir which has an K\ ll: of more 
than 500 /(M compared to 1 /<M of thymidine. 1 '' 

Kquation 1 shows that, for this series of compounds, 
affinity for the enzyme is affected by both the hydro 
phobicity of the 5-substituted side chain and its intorac 
lion energy with the enzyme. According to this model, 
strong affinity results from the binding of hydrophobic 
compounds which are also characterized by their ability 
lo interact strongly with the enzyme. Therefore, from 
the viewpoint of structure-based drug design, it may 
be of interest, to investigate which geometric or eiec 
ironic factors do influence either the interaction ener- 
gies of the 5-substitut.ed subst rates or their hydropho- 
bicities. 

(A) Factors Influencing the Interaction Energy, 
(i) Steric Complementarity with the Active Site 
Pocket. A general observation is that the intei action 
energy between substrate and enzyme decreases as the 
size of the ihienyl of furanyl substituent increases 
(Table 4. distance dependent dielectric with 3 A residue 
cutoff). Comparing, for example, the 2-thieny! ana- 
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logues, ii is clear that the smaller, unsubstituied 
heterocyclic thiophene ring 1 His the active-site pocket 
belter than ihe larger, bromo-substitutod thicnyl com- 
■jnunds Z 4. The same can be said For the substrates 
in the / S'uranvl series (10 > 12 > 11). or fur the 
isuxazuM scries ol compounds (14 > 15). Therefore, 
j substituted hcteiocwiir ring systems with small or 
no subst it nl ions aic preferred in terms of favorable* 
interact ion energv. 

(ii) HIeetrostatic Complementarity with the Ac- 
tive-Site Pocket, l'he interaction energy of isoxazolyl 
14 or funnel derivative 10 is more favorable than the 
corresponding value of thiophene compound 1. Figure? 
9 compares \he predicted binding conformation of 
compound 14 with that of 1. From this figure it can be 
seen thai in the case of the isoxazolo ring 14 a stabiliz- 
ing interact inn is formed due to the location of the two 
electronegative atoms 01 and NX in a region of positive 
electrostatic potential In the case of the thiophene ring 
i this stabilizing interaction is clear! v absent because 
c,i tin- lock of anv elect i onegai ive atoms in this region 
of positive potential, i'hus. the introduction of elec- 
tronegative atoms in the hcieroi \clic ring system at a 
position that allows them to interact with the positive 
potential at the back of the pocket may stabilize* the 
into! action between substrate and enzyme to a certain 
extent 

(B) The Counterbalance between Hydrophobic- 
ity and Interaction Energy. In t he case of Substi- 
tuted X'-deoxvuridine nucleosides, improved binding 
affinity for the* enzyme can be obtained by ameliorating 
the compound's hydrophobieity by introducing halogen 
or methyl groups. However, the introduction of large 
groups onto the heterocyclic system often leads to a 
decrease in interaction energy with the enzyme. This 
in turn translates into decreased affinity so that the 
final affinity will remain somewhat unaffected. For 
example. Hgurc \\ shows there is no remarkable differ- 
ent e in mcasure( ! aflin.u between the thiophene sub- 
sume 1 and lis '/ methvl. 2 chlorine, or Xbiomine 
analogues (b. 5. >: 2 f espect i\ eh ) . I lowover. inl.roduc 
i iun ni t i \ Murine 1 Z m hi inniiu* ;iicmi 11 onto a fui anyl 
ring 10 leads to a decrease in alliniu of at least 1 order 
of magnitude. I his is because the gain in affinity due 
to improved hvdi ophobicit \ of 1 1 and \Z is counlerbal 
a need to n much target extent by I he loss of affinity due 
io the significant decrease in interaction energy as a 
consequence of sum ic clashes between the side chain of 
Arg I IK* and the halogen subst intent. 



Arg-163 




Gln-125 

Figure 9. The electrostatic potential at the solvent-accessible 
surface of the active-site pocket. Red is positive potential: blue 
is negative. The C atoms of the isoxazule compound 19 are 
coloied jjray. and the thiophene analogue carbons 1 are coloied 
cyan. 

A similar conclusion can be drawn about the iut.ro 
duct ion of electronegative atoms in the five mombered 
ring systems. Although the replacement of ring carbons 
with electronegative atoms often improves the interac- 
tion energy due to the favorable interaction with the 
positive electrostatic potential near the back of the 
pocket, the resulting structures are characterized by a 
decreased hydrophobieity. This is exemplified by com- 
paring the affinities of the 2-furanyl substrate 10 with 
those the isoxazolyl compound 14. The small gain in 
affinity due to the increased interaction energy of 14 
compared to that of 10. is completely counterbalanced 
by the significant loss of hydrophobieity of 14 compared 
to that of 10. 

It is clear that structural changes aimed at improving 
the 5 substituent's hydrophobieity are often rounterbal 
anted bv decreased interaction energies and vice vorsn. 
Therefore, for the design of improved o -substituted 2' 
deoxyuridine analogues, it is the ding designer's dial 
lenge to search for that molecule which fortius the best 
compromise between hydrophobirit y and structural 
complementarity to the enzyme. 

Validating the Model. In order to assess the 
significance of the above derived correlation between 
interaction energy, hydrophobieity. and affinity (eq I), 
its success in predicting the ICw of 13 (pyrrol 1 yl) 2' 
deoxyuridine (18) (Figure 10) was investigated. This 
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Figure 10. I he strui lure ol a (pvi idI lyl)-2'-dooxyuridine (18) 
used lo validate U u ■ derived nn relation between interaction 
ee.eigv. hvurophuhicsiv and aflinity. 

compound had already heeu synthesized before, but its 
affinity for 1 ISV 1 TK had never been measured so far. 
:> (|»\ mil 1 \l)-?' deux\ uridine (18) was considered to be 
a suitable candidate because ol iis structural resem- 
blance v\i : . h t ujuptanuls like- those of Figure 3. A total 
suivatiori fire energy f«r 1 nieihylpyrrole of 9.1 kj/ 
nuil and. after docking compound 18 in the; active site 
of I ISV- 1 !'K. an into action energy of 304 k.J/mol were 
calculated. Following eq I. this corresponds to a 
predicted 11.' of i> «V1. which is in excellent agreement 
va ill) a measured value of 7 iiM. 

Conclusion 

riiymidine kinase from I ISV- 1 is a key enzyme in the 
metabolic activation of antiviral nucleosides. High 
affinity of these compounds for I he enzyme is required 
for efficient phosphorylation. Recently, a new series of 
:i substituted ?.' denxy uridine nucleoside substrates as 
ami 11S\ t agents lias been discovered, and many of 
thrse compounds show high affinity for the viral TK 
eu/vme. The availabiiit v of t hese data in combination 
v\ ith the : erenth solved structure of I ISV - 1 TK was an 
amartiw opportune. \ to investigate the structural 
fiiciors v\hic!i would influence the affinity of these 
compounds loi i he cn/\ me 

Modeling calculations indicated that the binding 
energetics of thymidine and 5-subst ituted 2'-deoxyuri- 
dine analogues are very similar. The majority of the 
binding energy arises from hydrogen-bonding interac- 
tions involving sugar and base moieties. Stacking 
interactions of the m m il base contribute the remaining 
binding forces. 

Calcuhn inns also indicated ihut the observed affinity 
differences ol the o substituted nucleosides are fine 
tuned by either I he hvdrophobicit \ of the 5-sub.st ituted 
side chain or its interaction energy with the enzyme. 
! Ivdi nphobii it v can be decreased b\ the introduct ion of 
halogen vubst itun nt^ whereas the interaction energy 
i-iin in- implied bs the introduction of elect conceal ive 
atoms a: uell defined positions, m by keeping the 
f>eometrv of ! he li\ e meinbri rd hewi orycle small. 1 low 
e\er. foi ibis sei ies ol i) substituted nucleosides, struc- 
tural modifications aimed at improving the compound's 
hydrophobic itv are often counterbalanced by decreased 
interaction energies, and vice versa, which makes the 
design of potent a suhsi ituted substrates an interesting 
challenge The present study also explains why 5-(bro- 
moviuN 1) ^'-deow uridine is one oi the strongest binding 
compounds, and it. finally, gives a qualitative explana- 
tion for the first metabolic step in the activation of these 
class of compounds. 

Fxperimental Section 

The ItSV i TK Crystal Structure. I'he crystal structure 
o! 1 ISV i i K in (i.!ui)ir\ w\ih ihviiucliiu- was used as starling 



structure." All ( rvstallugraphie waters were lomoved. liH 
PLUS- was used to provide an analysis of the hydrogen 
bonding pattern around the Asn. Gin. and His side chains and 
if necessary to exchange their N/O atoms where the hydrogen 
bonding pattern suggested it. Hydrogen atoms wore added 
at geometrically reasonable positions using MacroModel/" 
Polar hydrogens oi' Ser. Tyr. and Thr residues were reposi 
tioned manually based on geometrical and hydrogen-bonding 
criteria. Ail His residues were assumed to be neutral. Geo- 
metrical and hydrogen bonding criteria were used to decide 
between the two tautomeric configurations. Three Mis resi 
dues were modeled with their Nt atom protonated (58. M2. 
213). and the other five His residues were protonated ai No 
(105, 164. 283. 323. and 351). 

Charge and Force Field Parameter Development. The 
AMBER* force field 24 was used. Partial atomic charges for 
the substrates were calculated using the standard piece 
dure of fitting the electrostatic potential at the STO-3G 
level of theory. With the exception of torsion angle and 
some bond and angle parameters, all necessary parameters 
fur the substrates were provided by the AV1BKK' lone 
field. Missing bond and angle parameters were 1 derived from 
comparison with the crystal structures, while iheii lone 
constants were taken from comparable bonds or angles from 
the AMBFK* force field. In cases where no torsion parameters 
wore available, appropriate parameters were derived by 
comparison with torsion profiles calculated at the 6 31G* level 
of theory. However, because we were interested not in 
molecular dynamics but only in energy minimization calcula 
tion.s. little or no attempts were made to accurately model the 
height of barrier between local energy minima. A modified 
AMBFR* force field file in MacroModel fori rial containing 
partial charges and modified force field can be obtained from 
the authors upon request. 

Hnergy Calculations. Kncrgy calculations on enzyme 
substrate complexes were performed using the Batch VI in 
module of MacroModel 5.0."' No explicit waters or other 
approaches to model solvent were used because the substrates 
in die active site pocket are totally shielded from bulk solvent. 
Furthermore, nearly all solvent -accessible residues were kepi 
constrained to their crystal coordinates (see further), so that 
possible conformational changes due to the absence of solvent 
during the calculations were avoided. Aflei the substrates 
were docked in the active site, the whole complex was energy 
minimized until the energy gradient dropped below 0.01 k.F 
mo! -A. Other minimization protocols, such as simulated 
annealing procedures by which the systems were slowly cooled 
from 300 to 10 K over a 300 ps time scale, followed by energy 
minimization, did not improve the results to any extent {data 
not shown). Similarly, medium temperature molecular dy- 
namics in conjunction with energy minimization to explore 
conformational space did not yield new binding conforrna 
lions.™ Kit her a distance-dependent or a constant dielectric 
was used to calculate electrostatic interactions using a non 
bonded cutoff distance of 10 A. Only part oi the enzyme was 
allowed to relax, and residues furthei than a specified disunite 
10. 2. 3. T (i. or 8 A) from tin - substrate were kept constrained 
it) their crysiallographic positions using a harmonic force 
constant of 1(100 kJ/mohA' 

Calculation of Solvation Free Knei^ies. i'he electro 
static component of the total solvation tree energy was 
calculated by solving die Boissou Bollzmann equation"' using 
a finite difference 'method as implemented by the Dell'hi 
program." 0 n ™ The solvent-accessible surface was calculated 
with the DMS program of MidasBlus -' using a probe radius 
of \A and 0 A for the solvent calculations [f 80) and vacuum 
[( • 2) calculations, respectively. Partial atomic charges were 
obtained from 6-31G** calculations using the Gi IKKPG pro 
cedure of GAMKSS.*' The reduced radii set of Sitkoff and co- 
workers was used for' the atomic radii to calculate the solvent- 
accessible surface. 20 This procedure should guaranty the 
calculation of solvation energies with an average accuracy 
better than 3 U/inol. 20 Other charge sets, such as the STO- 
3G-derived atomic partial charges of the AMBKR* force field, 
were also Hied but yielded much less satisfying correlations 
(data not shown) 
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llu- no5»|)o];ii rumponrni i>i (ho tola) solvation froc energy 
was fSiiniHicri from the sulwni Miressihlc .surface area using 
Manikin! Pauling radii. A linear least squares relationship 
between ihe nonpoiar component or the total solvation free 
eni»rg\ anil the sulveni accessible surface area was calculated 
hum 'experimental nlkaiic transfer free energies." 1 ' and this 
Irasi squares illation served to extrapolate the nonpoiar 
solvation energies "1 1 1 if loinpounds under sludv. 
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